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Abstract 

We discuss the averaging hypothesis tacitly assumed in standard cosmology. Our ap- 
proach is implemented in a "3+1" formalism and invokes the coarse graining arguments, 
provided and supported by the real-space Renormalization Group (RG) methods, in par- 
allel with lattice models of Statistical Mechanics. Block variables are introduced and the 
recursion relations written down explicitly enabling us to characterize the corresponding RG 
flow. To leading order, the RG flow is provided by the Ricci-Hamilton equations studied 
in connection with the geometry of three-manifolds. The possible relevance of the Ricci- 
Hamilton flow in implementing the averaging in cosmology has been previously advocated, 
but the physical motivations behind this suggestion were not clear. The RG interpreta- 
tion provides us with such physical motivations. The properties of the Ricci-Hamilton flow 
make it possible to study a critical behaviour of cosmological models. This criticality is 
discussed and it is argued that it may be related to the formation of sheet-like structures in 
the universe. We provide an explicit expression for the renormalized Hubble constant and 
for the scale dependence of the matter distribution. It is shown that the Hubble constant 
is affected by non-trivial scale dependent shear terms, while the spatial anisotropy of the 
metric influences significantly the scale-dependence of the matter distribution. 

SISSA Ref. 8/95/FM (Jan 95), (Improved Version) 
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1 Introduction 



A successful physical theory usually enables us to isolate some limited range of length scales, 
or select a not too big set of variables, to render the problem tractable and at the same time 
preserve its essence. Fortunately, in many circumstances it is not necessary to resolve the details 
associated with each scale, since generally phenomena of each size can be treated independently. 
For example, in hydrodynamics there is no need to specify the motion of each water molecule 
and yet waves can be described as a disturbance of a continuous fluid. 

There are also problems that have many various scales of length, that is to say, phenomena 
or processes where each length-scale's contribution is of equal importance. To handle them one 
has to take into account the entire spectrum of length scales, dealing with fluctuations of almost 
any wavelength and consequently many coupled degrees of freedom. For example, critical 
phenomena, turbulent flow, the internal structure of elementary particles and confinement in 



QCD belong to the above class of problems |50|. 



We argue here that also in gravitational physics we encounter a problem of a similar nature, 
namely, the so-called averaging problem in relativistic cosmology. 

Relativistic cosmologies rely on some form of Cosmological Principle. The latter is usually 
a smoothing-out hypothesis imposed a priori on the distribution of matter in the universe. A 
well known example is provided by the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 
which is assumed to describe the real universe. 

There are at least two reasons why one should be very careful in doing so. 

Firstly, observations of the distribution of matter have shown that the scale at which the 
background homogeneity is reached,Q is probably of the order of hundreds of Mpc, meaning 
that locally, i.e. at least up to this scale, the universe is quite lumpy. Consequently, the 
local geometry is very complex and its nature not very illuminating from the point of view 
of cosmology. Seen this way, the approach usually taken is to average out all the matter, i.e. 
to redistribute them in the form of a homogeneous perfect fluid and use continuous functions 
(e.g. matter density, pressure) in modeling the universe, assuming that they represent "volume 
averages" of the corresponding fine scale quantities. In doing so one tacitly assumes that such a 
smoothed-out universe and the real locally inhomogeneous one, behave identically under their 
own gravitation. However, it has been stressed in [14] that the above assumption, though 



usually taken for granted, is by no means justified. First of all, Einstein's equations are highly 
non-linear which is why any averaging process on them|^ is far from trivial in general. 

Secondly, and this is a remark of a philosophical flavour, most of the observational data 
are theory dependent, i.e. their meaning can be interpreted only by assuming a particular 
theoretical explanation. 

Thus it is of great importance that the very foundations of a cosmological model should be 
as sound as possible, in particular, they should be free from assumptions which may not be 
warranted. 



If this scale has been reached at all is still a matter of controversy. 
2 Moreover, additional care is required since a volume average of a tensor is not a covariant quantity, unlike 
scalars. 
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2 Coarse-graining in Cosmology 



A possible solution to the averaging problem would be to explicitly construct a procedure 
for carrying out the smoothing process in the full theory. Almost all existing attempts were 
concerned with the linearized theory, with a possible exception of |]] (see also p4| ), for a review 
0. 



sec 



In H a covariant smoothing-out procedure was put forward for the space-times associated 
with gravitational configurations which may be considered near to the standard ones, gener- 
ating closed FLRW universes. The procedure makes use of Hamilton's theorem about smooth 
deformations of three-metrics and is adapted for smoothing-out an initial data set for cosmo- 
logical solutions to the Einstein equations. While interesting in its own right, this approach 
seemed rather ad hoc and not yet capable of resolving the issues of actual limits of validity of 
the FLRW models in cosmology. 

Now, our hope is that there is a smart and simpler way to the heart of the problem, bor- 
rowing from the known theories and methods of statistical mechanics, based on the real-space 



Renormalization Group (RG) approach to study critical phenomena in lattice models [27, 48]. 
Although the usual renormalization transformations, invoking averaging over square blocks, are 
designed mostly having ferromagnetic systems in mind, there are many more problems suitable 
for RG methods. These are difficult problems where the reductionist approach fails and where 
the effective degrees of freedom of a physical system are scale dependent. The difficulty of this 
kind of issues can be traced to a multiplicity of scales and, moreover, there can presumably be a 
gross mismatch between the largest and smallest scales in the problem. The averaging problem 
in cosmology can be looked at and studied as belonging to precisely this kind of problems. This 
is our main objective in this paper. 

Some form of RG is active on any system where there are fluctuations present (they by 
no means do need to be quantum). This is so, since one can integrate the fluctuations out of 
the physical quantities of interest, e.g. the partition function, and depending on the "scale" 
up to which one is integrating the same emerging quantities are different. The functional rela- 
tions between them provide recursion relations between the physical parameters, the coupling 
constants, which characterize the physics at each scale, and this is precisely what RG is all 
about. 

Often a major step consists in finding a way of looking at things. Therefore we stress that 
the problem we face with the averaged description in cosmology is effectively a question of how 
a system behaves under changes of "scales" . As such it is most naturally addressed using RG 
approach, understood here rather as a general strategy to handle problems of multiple length 
scales^] enabling us to extract the long distance behaviour of the system by making the scale 
successively coarser. In cosmology, we have the curvature inhomogeneities and to consistently 
tackle this problem, we will have to consider a procedure operating on the metric, not only on 
or apart from the matter present. 

To provide even more support in favour of the presented above idea, that RG arguments 
might be applicable to the problem of the inhomogeneous universe, let us notice that one can 
also be guided by scaling ideas. Scaling is exhibited (approximately or exactly) by many natural 

3 Although the renormalization procedure might seem purely formal there are important physical ideas behind 
it, namely, that of scaling and universality. 
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phenomena and mathematical models. The universe, namely, the distribution of matter in it 



does exhibit certain scaling properties |7], 26], of which the power law behaviour of the two-point 
correlation function for galaxies, clusters and quasars, is a fair example. Scaling, on its own 
right, is deeply understood within the underlying mathematical scaffold which is RG. These 
are hints therefore that one can regard the universe as a gravitational dynamical system not far 
from criticality (understood intuitively by analogy with e.g. a ferromagnet). Later one can also 
try to qualify the precise nature of the critical behaviour within the phase transition context, 
but this will not be of our concern in the present paper (apart from a simplified example in 
section I 



In passing, let us note that interesting results were obtained in [|T^] on spherical scalar 
field collapse and in 0] on axisymmetric gravitational waves collapse, which show a surprising 
scaling and critical behaviour reminiscent of that found in many (second-order) phase transition 
phenomena in condensed matter. They could be described, we believe, using RG approach 
again, treating the Einstein equations as generating a RG flow on the space of initial data. 

The real space renormalization techniques are mostly applicable to discretized models, based 
on a lattice. Therefore we now turn to describing a suitably discretized manifold model we are 
going to work with. 

2.1 Discretized manifold model 

The approach taken is that of a (3+1) formulation of General Relativity (GR) ||]. Let us 
suppose we have a differentiable, compact three-manifold (without a boundary) A4. Generally, 
in this case we will always assume that these manifolds possess certain natural constraints on 
their diameter and a suitably defined notion of curvature. The point of this requirement is that 
the manifolds, or more precisely the riemannian structures, can then be classified according to 
how they can be covered by small metric balls (to be defined later). Moreover the space of 
riemannian structures has some remarkable compactness properties. This is a classical result 
obtained by M. Gromov ^C] ]. On a set of riemannian structures it is possible to introduce a 
distance function, the Gromov distance, which roughly speaking enables one to say something 
about how close particular manifolds are to each other. For the riemannian manifolds which 
can be considered close to each other (in the sense of Gromov distance) it is possible to cover 
them with the balls arranged in similar packing configurations |Tot| . 



In order to define such coverings [21], let us parameterize the geodesies by arc- length, and 
for any point p G M. let cIm(x,p) denotes a distance function of the generic point x from the 
chosen one p. Then for any given e > it is always possible to find an ordered set of points 
{pi, . . . ,pn} in M, so that pl| 



i) the open metric balls (the geodesic balls) B_M(pi, e) = {x £ M.\ dM{ x -,Pi) < e}> i = 1, ■ ■ ■ ,N, 

cover A4; in other words the collection {pi, . . . , p^} is an e-net in A4. 

ii) the open balls Bj^ipi, e/2), i = 1, . . . , N are disjoint, i.e., {p±, . . . ,pn} is a minimal e-net 

in M. 

It is fair to say that as a consequence of the compactness properties of the set of riemannian 
structures that we consider, for each "length scale e" there exists a finite number of "model" 
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geometries which describe, with an e— approximation, any given riemannian geometry. Namely, 
given a ball of a certain radius > e in any riemannian manifold (with suitable restrictions on 
their volume, diameter and sectional curvature, as we have said earlier) there exists a ball 
metrically similar (up to an e scale) in one of the model geometries which does not retain the 
details of the original manifold on scales smaller than e. Roughly speaking, e is a measure of 
the typical curvature inhomogeneity with respect to the model background. Let us stress that 
this is a highly non-trivial result, in the sense that the metrical properties of the manifolds 
from an infinite dimensional set are, up to an e scale, described by the metrical properties of 
just a finite number of model riemannian manifolds. 

The e— nets underlying the balls-coverings precisely provide the discretized manifold model. 
This coarse-graining of a manifold according to Gromov is the most natural coarse-graining 
one can think of, pertinent for manifolds with a lower bound to the sectional curvature. This 
assumption does not limit the generality of our analysis which is basically motivated by a 
concrete physical problem, whose nature allows us to deal from the beginning with manifolds 
that are already in a certain sense quasi-homogeneous (Cf. comments on the solvability of the 
Ricci-Hamilton flow, later on). 

In what follows, when speaking of balls we will always mean geodesic balls here. 

2.2 An empirical averaging procedure 

We assume that we have chosen a particular space-like hypersurface E of the four-dimensional 
manifold, on which the average of a scalar function / : E — > 1R is given as 

where vol(H,g) = J^dfig and (i g is the riemannian measure associated with the three-metric 
g of E. Since at this stage we simply wish to put forward a few elementary geometrical 
considerations, we do not specify yet the choice of the hypersurface E and we do not attribute 
any particular physical meaning to the function /. 

If the geometry of E is not known on a large scale, we cannot take ([l]) as an operational 
way of defining the average of /. From a more pragmatic point of view, supposing that we can 
only experience geometry in sufficiently small neighborhoods of a finite set of instantaneous 
observers, it makes much more sense to replace (|l]) with a suitable average based on the 
geometrical information available on the length scale of such observers. 

For simplicity, given a finite set of instantaneous observers, located at the points xi, . . . ,xn £ 
E, we may assume that these susceptible to observation regions are suitably small geodesic balls 
of radius e, scattered over the hypersurface E so as to cover it. In other words, we assume 
that {x\, . . . ,xn} is a minimal e— net in E. Further, we denote by U e the corresponding set of 
geodesic balls {B%(xi, e)}, i = 1, . . . , N. Then we can bound (|IJ) as 

Eft vol(B h ,e/2) Eh vol ( B h,e) 

< J >e/2 7777 \ < < J >E( 9 ) << J >e 777; 7 , (^) 

where 
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J2i vol{Bi,e) 

and 

, f . _ Di/(Bj, e /2) f d ^9 (A , 

7 J2iVol{Bi,e/2) 

and where we explicitly indicated the dependence of the averages on a particular "covering". 
This suggests to consider < / > e as a suitable scale dependent approximation to < / >. 

There are certain problems lurking that we have to clear up. Obviously, there are "un- 
wanted" details affecting the average function over the discretized manifold, as given by its 
partition with a collection of geodesic balls, the immediate one being the underlying discretiza- 
tion. The important question to ask is what happens to the average when we change the length 
scale. Depending on whether we are actually increasing or decreasing it, respectively less or 
more details of the underlying geometry, will be felt by the average values. The natural philos- 
ophy is that over scales big enough no details should be discerned since the homogeneity and 
isotropy prevails. This is the reason why on constant curvature spaces averaging is well defined 
since there one can move the balls freely and deform them, but by so doing no new geometric 
details that measure the inhomogeneities will be felt in the averaged values of quantities we 
are interested in. 

A natural question to ask now is then how the geometry, specifically curvature inhomo- 
geneities, should depend on scale so that the average over the balls is scale independent, or 
equivalently, how do we have to deform the geometry in order to achieve the scaling limit when 
size of the balls matters no more? 

First, we give some details on how we calculate the averages according to ([|). In order to 
do this, we employ a preferred system of coordinates on the set of balls {B{} given by the local 
diffeomorphism 

exp x : T x Ti — > S (5) 
i.e., we make use of the exponential mapping 

<Pi = exp | exp _x B i =£) i : Di^ Bi, (6) 

where Di = D(xi, e) is the ball in T^£.[^] 

On Di we use polar coordinates and pull-back the riemannian measure accordingly, namely, 

<p*{Hg) = 9(t,Xi)dt®dxi, (7) 



4 The transition from (£,<?), at a given moment of time t a , to a tangent space parallels the prescription in 
riemannian geometry for measuring []43|. Let U be a given neighborhood of an instantaneous observer in E and 
suppose it is so small that there is a neighborhood A of 6 T^E such that exp^. : A(c T x T,) — > C/(C E) is a 
diffeomorphism. One can then replace the considerations in (U, g) by those in A (with the riemannian measure 
pulled back) via exp" 1 . Namely, we can say that an instantaneous observer in (E,g) observes the universe with 
the help of the exponential mapping, which just means projecting structures /.from an open neighborhood U C E 
of x by exp" 1 and treating them as structures on T^E. 
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where dxi denotes the canonical measure (euclidean volume form) on the unit sphere D{xi, 1) = 
Sf C r Xi S and where dt is the Lebesgue measure on 1R (t > 0). 

For t small enough one can prove Puiseux' formula 

9(t,x)=t n - 1 (l-±r(x)t 2 + 0(t 2 )), (8) 

where n = dim's and r(x) is the Ricci curvature Ric(g) (at the point x). 
Using this result we have 

/ fdfx g = [ f9(t,x)d Xi dt. (9) 

JB{xi,t) J Si 

Let us consider the asymptotic expansion with respect to t 

/ fdn 9 = u n t n [f( Xl ) + —^--(A/^) - ^/(s,)) + 0(i 2 )], (10) 
JB(xi,t) 2(n + 2) 3 

where uj n is the volume of the unit ball of IR n , R is the scalar curvature at the center of the 
ball and A the Laplacian operator relative to the manifold. 

Substituting / = 1 in the above formula we get the asymptotic expansion of the volume of 
a geodesic ball: 

vol(B( Xi ,t)) = u n t n (l - T^r^t 2 + 0(t 2 )). (11) 

o(n + 2) 

These standard formulae are what we need in order to calculate how the average value 
behaves when we change the scale. 

Since we are interested in discussing how < / > e behaves upon changing the radius of the 
balls {B(xi, e)}, let us consider the average < / > eo+ri , with rj a positive number with r]/e -C 1. 
According to the formulae recalled above, we can write 

* ; ^ s.[i-^S^i ' + 0({ '-° + " (12) 

where we have introduced somewhat simplified but otherwise obvious notation. Upon expand- 
ing this expression in 77, we get to leading order (C(e„) in e D , and 0([^] 2 ) in ^), 

< / >e D +„-< / >e + ~ ^< A/ > eo € 2 - + - 1 [< R > €o < f>e ~< Rf > ' e a W ~ , 

(13) 

where < / > £o is the average of the function / over the set of N instantaneous observers U to , 
(with similar expressions for < R > €o , < Rf > eo and < A/ > £o ). Thus, under a change of the 
cutoff we can write 



Zo^< f >, o+v \ V /e o =0 = ^~ 2 < A/ > eo + ^^y[< R >.< / >e a " < Rf >ej (14) 

to leading order. In the next sections we will discuss the consequences of these formulae and 
the connection of our averaging procedure with the Ricci-Hamilton flow. 
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3 The Renormalization Group view 



3.1 Block variables and recursion relations 

The real space RG technique is based on the recursive introduction of block variables. A method 
of "blocking" is trivial to introduce on regular lattices. In particular, in the case of the Ising 
model (2— dim) it consists of a subdivision of the spin system into cells, which supposedly 
interact in a similar way as the original spins. This can be done by introducing the block-spin 
variables via the majority rule or decimation procedure H. The behaviour of blocked lattice 
on large scales is equivalent to the behaviour of the original lattice corresponding to a different 
temperature. The slopes of the parameters' surface close to the critical fixed point determine 
the macroscopic characteristics of the model (see e.g. |]]). 

In the context of the problem we are considering, the application of RG method forces 
us to invent some analogue of Kadanoff 's blocking (block-spin transformation) applied to the 
geometry itself. This appears to be a difficult problem since in a general case when the geometry 
is curved the "lattice" itself takes on a dynamical role. Moreover, since the manifold £ we are 
dealing with is compact, (closed, without boundary), we must implement a Renormalization 
Group strategy in a finite geometry, and thus the relevant phenomena are here related to finite 
size scaling. Roughly speaking, the size of our manifold is characterized by a length scale, say 
L, (actually the volume, the diameter and possibly bounds on curvatures), which is large in 
terms of the microscopic scale (the radius of the typical geodesic ball coverings we shall use in 
the blocking procedure). A continuous theory, describing the (universal) properties of the field 
/ on S, arises when the correlation length associated with the distribution of / is large (of the 
order of L, or bigger). This being the underlying rationale, let us proceed and be guided by 
the formula (|l4|) . 

Let us consider our system as nearly infinite, i.e. the manifold £ divided by the collection 
of the geodesic balls of radius e m = (m + l)e D for m = 0, 1, 2, . . ., and where e Q , the chosen 
cutoff, is much smaller than the typical length scale associated with S, (this length scale can 
be identified with the injectivity radius of the manifold). Each ball will be labelled by k in the 
sequel and the original geodesic packing-covering is for m = 0. We can introduce this way a 
convenient notation for the integral of the generic function / over the ball B(xk, e m ) as 



which can be seen as the block variable since it allows us to eliminate ^from the distribution of 
the field / all fluctuations on scales smaller than the cutoff distance e m . We wish to emphasize 
that if the geometry of the ball B(xk,€ m ) is not flat, then the definition of ip m (k; /) can be 
interpreted as that of a weighted sum over a flat ball, namely, 



where the weight 9(t,x) is provided by Puiseux' formula (||). 

If we consider the covering of B(xk,e m+ i) induced by the geodesic balls {B(xj,e m )}, i.e., 




(15) 




(16) 
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the collection of N^irn^m + 1) open sets {B(xk, e m+ i)} n {B(xj,e m )}, then the above def- 
inition of block variables can be written recursively in terms of the values the function / 
takes correspondingly to balls of larger and larger radii. To this end, let us consider a par- 
tition of unity {C/i}/i=i,...,Jv fc (o,m+i)> subordinated to the covering of the generic enlarged ball 
B(xk, e m +i), induced by the geodesic balls {B(xh, £<>)}• Namely, a set of smooth functions such 
that: < £h < 1, for each h; the support of each £/j is contained in the corresponding B(xh, € ); 
and J2hCh(p) = 1, for all p £ 5(a; fc ,e m+ i). 

Under such assumptions, the block variables ip m (k; f) can be written recursively as 

JB(x k ,e ) 
JV h (m,m+l) 

i+ift/) = E ^m(fc;/). (17) 

Indeed, we have 

M h 'J) = / fdfig= J2 / 

JB(y h ,ei) k JB(x k ,e ) 

= E Mbf), (18) 
fc 

and 



JVj(l,2) 

ip2(j;f) = / /rfM fl = E / $hfdHg 

JB(zi,e2) I JB(y h ,ei) 



B(z 3 ,e 2 ) h JB(y h ,e 1 ) 

Nj (1,2) N h (0,l) 

k JB(xk,e ) 



= E E /, tehfdUg 



N h (l,2) 

= E ( 19 ) 

where we have exploited the fact that the functions have t,k£,i(p) = 0, except for a finite 

number of indices (A;,i), and J2i £fc£i(p) = 1) f° r an P £ X. The above expressions readily 
generalize for every m. 

More explicitly, we can write these block variables as 
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ip m+ i{k; f) 

N k (m,m+1) 



E 



B(x k ,e m +i) 



A(fZ h )(h)-R(h)f(h)£ h /3 \ , 4 
2(n + 2) ^m^^nJ 



(20) 



Notice that in terms of the block variables ip m (k; /) we can rewrite the empirical averaging 
(1) as 

T.k {m) vol[B(x k ,e m )} 

where N(e m ) denotes the number of distinct B{xj t ,e m ) balls providing a minimal e m -covering 
of the manifold E. Thus, when m is sufficiently large, the variation in < / > £m (ip m ) under a 
block transformation ip m (k) f) — ► ip m+ i(h; f) is given (to leading order) by 



(22) 



1 



n + 2 



<A(/&) > e e„ 



1 



+ 



1 



m + 1 3(n + 2) 



[< i? > £m < /a >e m " < Rfth >eje r , 



m 



The above choice of block variables brings out the coupling between averaging a scalar field 
over a manifold and the presence of fluctuations in the curvature of the underlying geometry. 

In order to be more precise, let us assume that the variables f(k) are randomly distributed 
according to some probability law P ({/(&)}) (later on we shall come back to this point with 
a definite prescription). Upon blocking the system and thus renormalizing the variables f(k) 
by increasing the scale size, the probability distribution P({f(k)}) induces a corresponding 
probability distribution on the variables ip m , viz., P({ip m (k; /)}). 

From equation (|l7| ) it is clear that if the geometrical properties of any two balls, B(xi, tm) 
and B(xj, em) (with B(xi, e m )r\B(xj, e m ) = 0), are not correlated then the corresponding block 
variables, ip m (i; /) and ip m {j]f), are uncorrelated. Such length scale L = e m characterizes 
the correlation (or persistence) length of the manifold (£,<?). It is a measure of the typical 
linear dimension of the largest ball exhibiting a correlated spatial structures. This correlation 
length can be seen in close analogy with the usual correlation length (usually denoted by £) in 
condensed matter systems. It depends there upon the coupling constants in particular upon 
temperature, and diverges to infinity at the phase transition point. 

Since g plays here the role of a running coupling, or if you prefer, of "temperature", the 
existence of a finite correlation length corresponds to a rather "irregular" , crumpled geometry 
(as seen on scales of the order of L), or, equivalently, a high temperature phase of our system. 
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According to the central limit theorem it follows, for m large enough, (e m S> L), that 
the block variables ip m (k;f), being the sum of uncorrelated random variables, are normally 
distributed, (let us say around zero, for simplicity), with a variance 



E P {^ m (i-f)) = N^ m h, 



(23) 



where Ep{. . .) denotes the expectation according to the probability law P({ip m }), \ is related 
to the variance of the variables f(k), R(k), and N^"'" 1 ^ denotes the number of e D -balls in the 



m-ball B(xi, e m ) 



Thus, irrespective of the details of the local distribution of the random variables f(k) and 
R(k), we can write for the distribution of {ip m (k; f )}i over (£,<?) 



dP({^ m })=\{ d^ m (k;f){27TN^x) 1/2 exp[-^ 2 m (k;f)/2N^x 



(24) 



This shows that by rescaling the block variables i^ m (k; f) according to 



<t>m(k;f) = 



-1/2 



i>m(k;f), 



(25) 



we get new block variables with a finite variance as m — > oo, and for random metrics (E, g) we 
can write 



dP({4>m}) = II [# m (A:;/)(2^)" 1/2 exp [-4> 2 m (k;f)/2x 



(26) 



The above remarks, paradigmatic of the real space Renormalization Group philosophy, 
show that the definition of a sensible blocking procedure, in our geometrical setting, consists 
of a transformation increasing the scale size, realized by passing from the variables f(k) to the 
variables ip m (k;f), (namely, by taking the average over all values of / in a larger and larger 
ball), followed by a rescaling obtained by dividing ip m (k; f) by a suitable power of the number, 
N^ e °' m \ of elementary e D -balls contained in the e m -ball considered, (for random geometries this 
power is 1/2). 



Following standard usage, and in order to arrive at an interesting geometrical notion of 
blocking, we assume that for a generic metric g this rescaling follows by dividing ip m (k; f) by 



jy(e ,m) 



, where ui m will in general depend on m. Thus, the rescaled blocked variables of 



relevance are 



4>m(k]f) 



l ip m (k;f). 



(27) 



The value of uj m will be fixed by the requirement that, as m — > oo, and for some (critical) 
metric g cr it, (in general for an open set of such metrics), such normalized large scale block 
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variables have a limiting probability distribution with a finite variance. Namely, 

m hm o P m ({0 m (A ; ;/)}) = Poo({0oo(/)}), 

E Poo (<PL(f)) = l. (28) 



Notice that if (£,<?) is a nice manifold, e.g. a constant curvature simply connected three- 
manifold, we are obviously expecting that the corresponding u m is 



w m (£,0) = l. (29) 
In general, we can assume that there is a set of critical metrics and a corresponding w, 



such that the above requirements (28) are satisfied. Such critical metrics are not necessarily 
constant curvature metrics and the corresponding u is not necessarily 1, (one may conjecture 
that 1/2 < uj < 1, as happens in the Renormalization Group analysis of many magnetic 
systems). At this stage, this is only a tentative assumption in order to arrive at an interesting 
concept of geometric Renormalization Group in our setting. Later on, we shall see that such 
assumptions are justified by exhibiting examples of such non-trivial metrics. 



4 Averaging matter and geometry 

Until now, our discussion has addressed mainly geometrical issues and the function / entering 
the cut-off dependent averaging < f > t was not specified. Now we wish to apply the obtained 
results to the averaging of the matter sources, namely, the matter density p, the spatial stress 
tensor s ab , and the momentum density J a , entering in the phenomenological description of the 
matter energy-momentum tensor with respect to the instantaneous observers comoving with 
E, viz. 

Tab = PUaUb + J a Ub + JbU a + S a b, (30) 

where u is a unit, future directed normal to the slice E, (u a u a = — 1). 

In fact, since in the present epoch the universe is mainly matter dominated, i.e. the pressure 
can be safely neglected, it would in principle be justified to start our analysis with 

T ab = pu a u b + J a u b + J b u a , (31) 

for the stress-energy tensor. Even the terms involving the momentum density can be eliminated 
if we pick up a sensible slicing St (the comoving frame) provided that the cosmological matter 
fluid is an irrotational fluid in equilibrium. In any case, wishing to maintain the following 
discussion to a sufficient degree of generality, without particular restrictive assumptions on the 
matter sources, we assume in our analysis that the matter energy-momentum tensor has the 
perfect fluid form 

Tab = pu a u b + J a u b + J b u a + g a bP, (32) 
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with a pressure p which is a priori not-vanishing, and where g a b is the three-metric of S, namely, 
9ab = + ^a^6j ffj^ denoting the space-time metric^. 

As argued in the previous section, smoothing-out the matter sources as described by a 
set of instantaneous observers (represented by the three-dimensional hypersurface E), means 
eliminating from the distribution of such sources on E all fluctuations on scales smaller than the 
cutoff distance e, leaving an effective probability distribution of fluctuations for the remaining 
degrees of freedom. The underlying philosophy being that this effective distribution has the 
same properties as the original one at distances much larger than e (i.e. for fluctuations with 
wavelengths much larger than e). 

In order to implement this idea along the geometrical lines discussed in the last section, we 
first need to better specify what we mean by the assignment of a collection of instantaneous 
observers (endowed with clocks) on E. Since we are adopting a Hamiltonian point of view, 
such observers are specified by the assignment, on the (abstract) three-manifold E, of the lapse 
function a and the shift vector field a 1 . The former provides the local rate of the coordinate 
clocks of such observers, while the latter is the three-velocity vector of the observers with 
respect to the set of instantaneous observers at rest on E. 

The macroscopic variables of interest characterizing the matter sources are in the present 
framework, the matter density p, and the momentum density J. Actually, it would make 
more sense to consider the cosmological fluid phenomenologically, described by a pressure p, 
a baryon number density rib, energy density p, specific entropy s, and temperature T; such 
variables being related by 

dp = ribdh — ribTds (33) 

and 



h = P ~^, (34) 
rib 

where h is the specific enthalpy. However, for simplicity, we shall in the sequel consider a 
barotropic fluid. 

Thus, the field / characterizing the matter sources, as described by the instantaneous 
observers on E, is given by 

f = ap + a l Ji. (35) 

Notice that 2ap{p) is the hamiltonian density of the fluid in the Hamiltonian formulation of 
Taub's variational principle for relativistic perfect fluids p9|]; we also assume that the dominant 
energy condition holds. 

5 The fact that we are adopting a stress tensor with a non- vanishing energy flux can be traced back to the 
observation that even in an exact fc = 0orfc = +l FLRW universe there can be a non-zero particle flux tilted 
with respect to the frame of the fundamental observers, even though T a t has the perfect fluid form. This is 
possible if the macroscopic quantities (e.g. a particle flux) are derived from kinetic theory with anisotropic 
distribution function f(x l ,p a ) 
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The averaging of the matter sources along the lines described in the previous section would 
then require considering a finite set of instantaneous observers, {x±, . . . ,x^}, located on £ 
and setting a standard for the cutoff distance e Q over which the (experimental) distribution of 
matter sources, (i.e. the probability that the matter variables, p(i) and J(i), conform to a given 
distribution p(i)dp, J(i)dJ), is determined. Then we proceed with the blocking prescription for 
eliminating unwanted degrees of freedom and consider the behaviour as the averaging regions 
become larger and larger. 

We may decide to treat the riemannian geometry of £ as uncoupled with the matter sources 
if we are simply interested in smoothing-out the sources, or if we wish to consider fluctuations 
of the matter as essentially uncoupled to the fluctuations in geometry. Namely, the curvature 
fluctuations appearing in the definitions of the block variables, il) m {k, f) and m (A:;/), can be 
thought of, under appropriate circumstances, as independent random variables with a given 
distribution. 

In general, however, the sources are coupled to the gravitational field, and we ought to treat 
the full dynamical system, the cosmological fluid plus the geometry of £, in the procedure of 
blocking. Moreover, we should bear in mind that without taking explicitly into account the 
backreaction of the geometry one cannot really provide a reasonable averaging procedure for 
the sources. 

In order to do so, we consider for a given cutoff distance e D , the variables p(k) and J(k) 
associated with a minimal geodesic ball covering {B(xk,e a )}^ =l . The change of the cutoff is 
naturally realized by considering balls {B(xi, e m = (m + l)e Q )} with m = 0, 1, . . .. The block 
variables, ip (k; p,J) and ip m (k; p,J), are defined according to ( |17| ) with / given by (|35|), This 
transformation can be seen as thinning-out the degrees of freedom which is at the heart of any 
coarse- graining. 

The original cutoff e Q is chosen to set the scale over which General Relativity is experimen- 
tally verified and it can be taken as, say, the scale of planetary systems. We can thus safely 
assume that the Einstein field equations hold on that scale. It is however rather impossible to 
provide a mathematical model of the distribution of matter in the universe going down to such 
fine scales; besides this task would be impractical. What one does instead is to use continuous 
functions assuming that they represent appropriately "volume averages". The results of such 
an averaging in an inhomogeneous medium obviously depend on the scale. The point is that if 
the Einstein equations hold on the scale where they have been verified (here taken to be that 
of the planetary scale), then they do not seem to hold a priori on larger, cosmological, scales 
that require averaging. To see this, notice that the Einstein tensor G^ u , calculated from an 
"averaged" metric 3™, cannot be equal to the Einstein tensor Guv which was first calculated 
/.from the fine-scale metric g^ v and then averaged. This is so due to the non-commutativity 
of "averaging the metric" with calculating the Einstein tensor being strongly non-linear in the 
metric components. 

Below we shall assume a Hamiltonian point of view. Then the probability P({p(k), J(fc)}) 
that the matter variables, according to the records of the instantaneous observers {B(xk, e); a(xk), a l (xk)} 
in £, take on some particular set of values, {p(k),J(k)}, is given by the equation 

P({Mk; P, J)}l,...,iv) = \ exp[-H({p(k),3(k)})}, (36) 
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where Z is a normalization factor, and H({p(k), J(k)}) is the Hamiltonian associated with the 
matter variables p and J. Namely, 



exp - < ap + a h Jh > to 

W<»<^ ^ E,,i P [-<^W < 37 > 

where J2 P ,j ■ ■ ■ denotes summation over all possible configurations of the matter variables, p 
and J, that can be experienced by the instantaneous observers {B(xk, e); a(xk), a l (xk)} in X. 

The validity of General Relativity at this scale implies that this Hamiltonian is a part 
of TiADM({p(k), J(fc)}), the Arnowitt-Deser-Misner Hamiltonian, associated with the data of 
the three-geometry g a ^ of S, its conjugate momentum 7r a 5 and of the matter variables p and 
J, evaluated in correspondence with the {B(xk, e)}-approximation associated with the net of 
points {xfc}. Namely, 

HADM{{p(k),J(k)})\ e =< aH(g, vr, G, p) + a h H h (g, vr, G, J) > e ^ vol{B(x h e)), (38) 

i 

where 

H(g,7r,G,p) = (detig))- 1 / 2 ^^ - ^«) 2 ] - (det(g)) 1 / 2 R(g) + Svr^et^)) 1 /^^, (39) 

and 

7~La{g, ft, G, J) = -2tt^ + levr^^)) 1 /^^, (40) 

and where the momentum conjugate to the three- metric g a b is given in terms of the second 
fundamental form K a b of the embedding of £ in the resulting space-time, as 

7r ab = {det{g)) l ' 2 {K ab - K c c g ab ). (41) 



We wish to recall that the lapse and the shift appear in the Hamiltonian as arbitrary 
Lagrange multipliers (their evolution is not specified by the equations of motion), and as such 
they enforce the constraints 

H{g,*,G,p) = (detig))- 1 / 2 ^^ - i«) 2 ] - {detig)) 1 ' 2 R{g) + S^detig^Gp = (42) 

and 

H a (g, vr, G, J) = -2ir b a . b + IQ^detig^GJa = 0. (43) 



As is well known, these constraints are related to the invariance of the theory under the 
(four-dimensional) diffeomorphism group of the space-time resulting from the evolution of the 
initial data satisfying them. The momentum constraint, TC a (g,ir,G,J) = 0, generates the 
(spatial) diffeomorphisms into S, while the Hamiltonian constraint, 7i(g, tt, G, p) = 0, generates 
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the deformation of the manifold S in the resulting space-time, i.e. the dynamics. (Notice 
that such deformations can be interpreted as four-dimensional diffeomorphisms only after that 
space-time has been actually constructed). 

The Hamiltonian Ti.ADAi({p(k) , 3 (k)}) , apart from the lapse and the shift, depends on the 
three-metric g a b on £ and its conjugate momentum ir a b, the gravitational coupling G, and on 
the configuration which the matter variables, p(k) and 3(k), take on the set of instantaneous 
observers, {B(xk, e); ct(xk), d l {xk)}, chosen to describe the distribution of matter at the given 
length scale e. 

The basic question to understand is how the block transformations, 

{<f) m (k;p,3)}^{<f) m+ x(k;p,3)}, (44) 

followed by rescaling, affect the Hamiltonian associated with the matter variables p and J, and 
then discuss how this in turn affects the full Hamiltonian TCadm- 

Starting from the probability distribution P({ip (k; p, J)}i,...,jv) we can inductively define 
the probability distribution p( m+1 \{ip m+ i(k; p, J)}) of the block variables {ip m+ i(k; p, J)} (and 
of the corresponding rescaled variables {4> m +i (k; p, J)}). Since the block variables {ip m+ i (k; p, J)} 
are recursively obtained from the knowledge of the block variables at the m th stage, {ip m (k; p, J)}, 
the distribution p( m + 1 ) ({^ m+1 (k; p, J)}) only depends on the knowledge of P^ ({ip m (k; p, J)}). 
We can formally write 

p( m+1 \{^ m+1 (k;p,3)}) = Yl P (m) ({Mk;p,3)}), (45) 

{ipm(k;p,J)} 

where the sum is over the probabilities of all the configurations {ip m (k; p,3)} consistent with 
the configuration {ip m+ i(k; p, J)} of the block variables. 

As usual, this allows us to define the effective Hamiltonian for the matter variables after m 
iterations of the block transformation, according to 

P {m) ({Mk;p, J)}) = ^eM-H {m H{^ m (k;p,3)})}, (46) 

where 

Z im)= £ eM-H {m) ({^ m (k;p,3)})]. (47) 

{ip m (k;p,J)} 



Such H( m \{ifj m (k; p, J)}) are defined up to an additive constant term (e.g. ||). If we 
also stipulate, as is standard usage in the Renormalization Group approach, that the effective 
Hamiltonian H^ m+l \{^ m+l (k;p, J)}) takes the same functional form as H^ 71 ' ({tp m (k; p, J)}), 
i.e., if 



({tp m+1 (k; p, J)}i,...,jv) 



exp 




< a p<tn+l) + a h j(m+l) 






E P , j exp 


- < ap(™+V + a h J { ™~ 


hl) > 1 



(48) 
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then this indeterminacy can be transferred to the effective matter variables p(' m+1 ) and j( m + 1 ). 
in terms of which H^ m+1 ^ ({ip m+ i(k; p, J)}) is defined. 



It is immediately checked that such effective matter variables are defined by ( |46|) up to the 
transformations 

p (m+l) _^ p (m+l) + ^ 

j(m+l) ^ j(m+l) +V) (49) 

where h and v are, respectively, a scalar function (sufficiently regular) and a vector field defined 
on E. 

Among the various normalization conditions that we may adopt, in order to avoid the 



indeterminacy connected with (49), the natural one comes about by requiring the divergence 
and the Hamiltonian constraints to hold at each stage of the renormalization. This being the 
case, upon coupling matter to the geometry, the lapse a and the shift a 1 maintain their role of 
Lagrange multipliers enforcing the (four-dimensional) diffeomorphism invariance of the theory. 

In this way, the renormalization is set up so that the constraints, if initially satisfied, will re- 
main satisfied at every renormalization step. The three-metric g, (and the second fundamental 
form), in its role of running coupling describing the interaction between matter and geome- 
try, is then governed through (Coulomb like) effects that are expressed by these constraint 
equations. In this sense, the Coulomb like part of the gravitational interaction is driving the 
renormalization mechanism, and such requirements imply that the full effective Hamiltonian, 
(matter plus geometry), takes on the standard ADM form pertaining to gravity interacting 
with a barotropic fluid at every stage of the renormalization. 

According to the remarks above, we can consider as an independent parameter in the 
(effective) fluid Hamiltonian ({ip m (k; p, J)}) the three-metric of the three-manifold 
E, whereas the matter density p^ and the current density j( m ) are at each stage connected 
to g^ and by the Hamiltonian and divergence constraints that hold at each stage. 

Then the effect of the renormalization induced by the blocking procedure {ip m (k;p, J)} — > 
{ip m +i(k; P, J)} and the corresponding rescaling, can be symbolized as a non-linear operation 
acting on the metric (g a ™ ) so as to produce the metric (g a ^ +l ^), i.e. 

(9T 1) ) = n 9a f), (so) 

whereas the renormalization of the second fundamental form K ao is generated by the lineariza- 
tion of (|50|), as will be discussed in section 5. 

This way, this renormalization transformation 1Z defines a trajectory in the space of rieman- 
nian metrics of E. One moves on such trajectory by discrete jumps. However, in what follows 
we shall replace such discrete dynamical system by a smooth dynamical system, describing 
renormalizations of the parameters (gab) and (Kah)- 

The deformation of the initial data set for the Einstein field equations, symbolically denoted 
by 1Z above in (|50|), realizes in fact a formal, (at least at this stage), mapping between the initial 
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data sets for the field equations. As we have seen above, this renormalization acts in such a 
way that at each its step the deformed data satisfy the constraints. The time evolution, in 
turn, of any such data set generates a one-parameter family of solutions to the field equations 
which interpolates between the initial inhomogeneous space-time to be averaged-out and its 
renormalized (deformed) counterpart. 

4.1 The Ricci— Hamilton flow 

In order to replace the discrete operator 1Z, describing the effect of renormalization of g a b, 
with a continuous flow we can start by discussing some geometrical implications of equation 



(14). According to the Renormalization Group analysis of the previous section, they follow by 
considering the average < / > t as a functional of the metric and thinking of the metric g a b as 
a running coupling constant, depending on the cut-off. In this connection, it can be verified 



that we can equivalently interpret (14) as obtained by considering the variation of < / > e 
under a suitable smooth deformation of the background metric g a b, rather than by deforming 
the (euclidean) radius of the balls {B(xi,e}}. 

As a matter of fact, we can equivalently rewrite the second term on the right hand side of 
(0) as 

<R> e <f> € -<Rf> e =-D<f> e -^, (51) 

where D < f > t -dgab/drj denotes the formal linearization of the functional < / > e in the 
direction of the symmetric 2— tensor dg a b/drj, and where 

(> ' U( ,,] ~ < R( V ) > gM - 2R ab (rj), (52) 



drj 3 

R a b(v) being the components of the Ricci tensor Ric{g(jj)), and < R(r]) > is the average scalar 
curvature given by 

< R( V ) >= f R( V )dfi v . (53) 

Indeed, the linearization of < / > e in the direction of the generic 2— tensor dg a b/drj is provided 
by 

D < f >e 'IT = \ < f9ab ]^ 9ab > '~\ <f >e< gab -^ 9ab >e ' (54) 

so (51) follows, given the expression (p2) for dg a b/drj. 



According to what has been said in the previous paragraphs, the effective distribution of 
matter sources, according to a set of instantaneous observers {B(xi,e)}, is characterized by 
the underlying three-geometry thought of as an effective parameter depending on the cutoff 
e. Since the Hamiltonian and the divergence constraints hold at each stage of the renormal- 
ization procedure, (they fix the effective Hamiltonians which otherwise are undetermined up 
to a constant factor), the renormalization of the matter fields is intrinsically tied with the 
renormalization of the three-metric. 
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The invariance of the long distance properties of the matter distribution, under simultaneous 
change of the cutoff e and the parameter g a }>, can be expressed as a differential equation for the 
effective Hamiltonian Ti.(p, J), (actually for the partition function associated with this effective 
Hamiltonian), namely, 

[-^ + /Ms)^]£exp[-W(p,J)] =0. (55) 

Recall that 7i(p,J) is explicitly provided by (in the given approximation) 

Hip, 3) =< ap{p) + a i J i > t J2 vol ( B (%h,e)), (56) 

h 

where the average < ap(p) + a l Ji > e is a functional of the three-metric g a b, here thought of as 
the running coupling constant. 

Thus, in order that the equation (|55|) is satisfied, it is sufficient that 



-e^r + PMa — ] < MP) + a l Ji > e = 0. (57) 
Ce og a b 



The Renormalization Group equation ( p7| ) states that increasing the cutoff length (i.e., the 
radius of the averaging balls) from e to e 7 e, while deforming the metric g a b by flowing along the 
beta function P a b(g) for a parameter-time 7, has no net effect on the long distance properties 
of the considered system. 

According to the above remarks, and equation (fll|), it can be verified that the beta function 
yielding for ([57]), defined by 

d 

e fe9ab = Pab(g), (58) 

is exactly provided by (p2f) , namely, 

Pa b (g) = = \< R(v) > Qabiv) - 2R ab (v), (59) 

where the parameter r\ is the logarithmic change of the cutoff length e. 



Since the manifold S is compact, (|52| ) has to be interpreted as a Renormalization Group 
equation in a finite geometry, and thus the relevant phenomena are here related to finite size 
scaling (see e.g. ||). A continuous theory, describing the (universal) properties of the cos- 
mological sources and of the corresponding geometry, may arise when the correlation length 
associated with the distribution of cosmological matter is of the order of the size of the under- 
lying manifold. 



The metric flow ( |52|) is known as the Ricci-Hamilton flow 1 22 ] , studied in connection with the 
quasi-parabolic flows on manifolds; quite independently it has been discussed in investigating 
the Renormalization Group flow for general cr-models (see e.g. [35] and references quoted 



19 



therein). The Ricci-Hamilton flow is always solvable [22| for sufficiently small rj and has a 
number of useful properties, (apart from being volume preserving which is simply a consequence 
of the normalization chosen), namely, any symmetries of g a b(i]o) are preserved along the gabi^l) 
flow for all 7] > r] , and the limiting metric (if attained) g ab = lim^oo g a b(r/) nas constant 
sectional curvature. Thus equation (|52j), with the initial condition g a b(0) = g a b, defines (when 
globally solvable) a smooth family of deformations of the initial three-manifold, deforming it 
into a three-space of constant curvature. 



4.2 Fixed points and basins of attraction 

Thus the point of the above discussion is that in order to arrive at a fixed point of the RG 
equation (|i"4|), the geometry has to be deformed according to the Ricci-Hamilton flow (52). 



In this setting of the problem, the Ricci-Hamilton equations appear naturally and in fact the 
approach proposed, enables us to attach a physical meaning to them within the coarse-graining 
picture. This element was lacking in Q where Hamilton's theorem appeared rather ad hoc. On 
the other hand, our approach demonstrates that the smoothing issue is deeply connected with 
the geometry and exhibits how this relationship works. First some general remarks. 

A fixed point is a point in the coupling constants space that satisfies 

9*ab=K(9*ab), (60) 

i.e. it is mapped onto itself by RG transformation. 

Under RG transformation length scales are reduced by a factor m + 1. Namely, for the 
block variables, the correlation length measured in units appropriate for them, L m (e m ), is 
smaller than the correlation length L Q of the initial system measured in units of e Q . The actual 
physical value of the correlation length L is of course unchanged by the process of blocking, 
thus L = L m (e m ) = L e , so 

L m = (61) 
m + 1 

Since L m < L Q , the system with Hamiltonian TiS m ' must be further ^from criticality than the 
original system. Thus, we conclude that the system is at a new effective reduced "temperature" 
(r? 1 . At a fixed point (|60[), L* can be only zero or infinity since then L* = L* j{m + 1). 

As is standard in RG analysis, we will refer to a fixed point with L = oo as a critical fixed 
point, when L = we will call it trivial. Each fixed point has its domain or basin of attraction, 
namely, the points in the coupling constants space in such a basin necessarily flow towards and 
end up at the fixed point, after an infinite number of iterations of RG transformation. 



Let us, for the purpose of clarity, employ for a moment the ferromagnetic analogy. In this 
case, for a system exhibiting a phase transition, there are two attractive fixed points. One 
is the high-temperature fixed point which attracts each point with T > T cr n in the coupling 
constants space, and it corresponds to the effective Hamiltonian for the system as T — > oo. In 
this phase the variables assume random values and are uncorrelated. Upon a sensible blocking 
of such a system the probability distribution of the block variables remains unchanged. 

The second fixed point is the low-temperature fixed point which is the effective Hamiltonian 
for the system when T — > 0. This corresponds to a system in a complete spins alinement and 
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the block variables are ordered then. Every point corresponding to T < T cr j t is eventually 
attracted to this fixed point. 



Across the Hamiltonian space there should exist then a surface, the so-called critical surface, 
which separates the effective Hamiltonians flowing to the high-temperature fixed point from 
those flowing to the low-temperature fixed point. Notice that the word surface here has rather 
a heuristic meaning. Indeed, the set of metrics, separating those flowing towards a "low- 
temperature" fixed points ^from those flowing towards a "high-temperature" fixed point, has 
quite a complex structure whose understanding is deeply connected with some, yet unsolved, 
conjectures in 3-D-manifold topology, (see below for more details). 

If we now choose to start with a point on the critical surface, then upon RG transformation 
it will stay within the critical surface. There is a possibility that, as the number of iterations 
of RG goes to infinity, the Hamiltonian will tend to a finite limit 7i*. This point is the critical 
fixed point and within the critical surface it is attractive, (this is roughly speaking the basic 
mechanism for universality), along the direction out of the critical surface it is repulsive. This 
fixed point is related to the singular critical behaviour of the system due to the fact that all 
points in its basin of attraction have infinite correlation length. The simplest case is when the 
fixed points are isolated points, but it is also possible to have lines or surfaces of fixed points. 

With these preliminary remarks along the way, let us discuss how some of the above general 
characteristic of the Renormalization Group flows find their proper counterpart in the particular 
situation we are analysing. 

In the previous paragraph we showed that it was possible to replace the discrete operation of 
increasing the scale size of the observational averaging region by a "transformation" , smoothly 
deforming the background metric g a b, which turned out to be the Ricci-Hamilton flow. In 
this setting, following the example above, we would like to adopt the fundamental hypothesis 
linking RG to the critical phenomena, namely, the existence of a "critical" metric (on the 
critical surface) and of a "fixed point" metric g* b , such that 



In ( |62| ) g* b is a mathematical object invariant under RG and we assume that g c £ h %t repre- 
sents the physics (in a sense to be clarified further) of curvature fluctuations of a manifold at 
its critical point ("critical manifold"). Since the Ricci-Hamilton flow can be interpreted as a 
dynamical system on a set of closed riemannian manifolds, we can adopt the following inter- 
pretation of ([62]) . Suppose, that we look at our manifold through a "microscope" and are able 
to discern the curvature fluctuations down to a size e m . Then T^ 77 ^ represents the operation of 
decreasing the magnification factor, by m say, i.e. the sample seen appears to shrink by this 
factor. We have to assume that the system is sufficiently large so that the edges of the sample 
will not appear in the view. The hypothesis (|6^) states that if we decrease the magnification 
by a sufficiently large amount, we shall not see any change if we decrease it even further. 

We are going to describe to what extent such hypothesis holds. 




(62) 
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We already said that the Ricci-Hamilton flow (|52|), while always solvable for sufficiently 
small 77 f?2l , may not yield for a non-singular solution as i] increases. Hamilton noticed that 
there are patterns in the kind of singularities that may develop. Typically, the curvature blows 
up but in a very regular way (e.g. for S 1 x S 2 with the standard symmetric metric). This has led 
him to a research program which, roughly speaking, amounts to saying that any three-manifold 
can be decomposed into pieces on which the Ricci-Hamilton flow is global and thereby, each of 
these pieces can be smoothly deformed into a locally homogeneous three-manifold. Singularities 
may develop in the regions connecting the smooth-able pieces, but such singularities should be 
of a finite number of types and all of a rather symmetric nature (namely, if they are blown up, 
they should be associated with symmetric manifolds as e.g. S 1 x S 2 [pi] , ^qj). 

It may be said that Hamilton's program is an analytic approach to prove Thurston's con- 
jecture, which claims that any closed three-manifold can be cut into pieces, such that each of 
them admits a locally homogeneous geometry. The rationale, underlying this analytical pro- 
gram towards Thurston's geometrization conjecture, lies in the above nice structural properties 
of the Ricci-Hamilton flow. Several steps are involved in this program. Let us briefly recall 



them since, even if some of them are yet unproven, they shed light on the assumption (62) and 
on (]52|) when interpreted as a Renormalization Group flow. 



The first step is the assignment of an arbitrary metric g on the three-manifold S. In the 
Renormalization Group approach, this corresponds to picking up a vastly inhomogeneous and 
anisotropic geometry for the physical space, (here equivalent to the high-temperature phase). 
Such a choice may not conform to the actual quasi-homogeneous three-geometry of the physical 
space as is experienced now. However this quasi-homogeneity, in our opinion, may be related 
to the possibility that the actual universe is near criticality, a circumstance that we want to 
discuss rather than assume from the outset. 

The second step is to deform this metric g via the Ricci-Hamilton flow (^). In general, 
this flow develops local singularities which should be related to the manifold decomposition in 
Thurston's conjecture. Away from each of the local singularities it is conjectured (but not yet 
proved) that the Ricci-Hamilton flow approaches that of a locally homogeneous geometry in 
each disconnected piece. 

This picture may be consistent with our blocking procedure yielding for (^) as the Renor- 
malization Group flow. According to the analysis, carried out in section 3.1, a highly inhomoge- 
neous and anisotropic geometry (S,<?) can be characterized by minimal geodesic ball coverings 
{B(xi, e )}, whose balls are to a good approximation largely uncorrelated when seen on a suit- 
able scale. This means that the values of the scalar curvatures, R(i), evaluated at the centers 
of the balls of the covering, are random variables. Upon enlarging the balls and rescaling, cer- 
tain regions of the manifold might be such that correlations develop among the corresponding 
R(i), whereas in other regions, no matter how we enlarge the balls and rescale, the R(i) will 
remain independently distributed random variables. The former regions, exhibiting a persis- 
tence length, are those that should approach a locally homogeneous geometry under blocking 
and rescaling (i.e., under the flow (0)). The latter regions should rather develop singularities 



under (52), since no matter how much we block and rescale, the curvature will maintain its 



white noise character. 
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The third step is to study the behavior of (|52[) for the locally homogeneous geometries, for 



this accounts for the structure of the critical set of (p2[). This has been accomplished [25], and 
one finds that, depending on the initial locally homogeneous geometry, the Ricci-Hamilton flows 
either, (i) converges to a constant curvature metric, (ii) asymptotically approaches (as 77 — ► oo) 
a flat degenerate geometry, of either two or one dimensions (pancake or cigar degeneracy) , with 
the curvature decaying at the rate 1/rj, or, (Hi) hits a curvature singularity in finite time, with 
this singularity being that of the Ricci-Hamilton flow for the standard metric on 5 2 x S 1 . Note 
that constant curvature geometries always occur whenever the manifold can support them, (in 
dimension three, constant curvature manifolds and Einstein manifolds are synonymous). It is 
also quite interesting to note that the Ricci-Hamilton flow of homogeneous metrics usually (with 
a few exceptions) tends to approach or converge to the maximally symmetric homogeneous 
metric in the class considered (see [25] for details). 



Interpreting ( |52| ) as the Renormalization Group flow, it follows that locally homogeneous 
geometries evolving under ( |52"|) towards an isolated constant curvature manifold, are sinks 
describing a stable phase of the corresponding cosmological model. For instance a FLRW 
model (with closed spatial sections) characterizes such a phase. Non-isolated constant curvature 
manifolds (e.g. flat three-tori) provide less trivial examples of limiting behaviour of (|52"|), (see 



e.g. [23, [ll|]). The locally homogeneous manifolds non-admitting Einstein manifolds (i.e. there 
are no left-invariant Einstein metrics on the group SL(2, M)), provide even more interesting 
behavior. In this case, a metric renormalized under the action of (^) develops degeneracies 
and one gets, in our setting, an effective cosmological model with spatial sections of lower 
dimensionality. 



All these limit points of (|52|), either fixed or not, have their own basins of attraction. Rig- 
orously speaking, they are the sets of three-metrics flowing under ([52]) to the respective limit 
points, briefly discussed above. There are nine such basins of attraction, corresponding to the 
nine classes of homogeneous geometries that can be used to model (by passing to the univer- 
sal cover) the local inhomogeneous geometries on closed three-manifolds. By labelling these 
classes according to the minimal isometry group of the geometries considered we distinguish 
the following basins (here we follow closely the exposition in [|25|]): 

(i): The 1R 3 -Basin. It contains all three-metrics flowing towards the homogeneous flat 
]R 3 metrics. This basin is eventually attracted by flat space, (flat tori, when reverting to the 
original manifold rather than to its universal cover). 



(ii): The SU (2) -Basin. It contains all three- metrics flowing towards the three-parameter 
family of homogeneous SU(2) metrics. This class admits Einstein metrics, in particular the 
round metrics on the three-sphere. This basin is exponentially attracted to the round three- 
sphere, (modulo identifications). It is the basin of attraction yielding for closed FLRW cosmo- 
logical models. 

(Hi): The SL(2, IR) -Basin. It contains all three-metrics flowing towards the three-parameter 
family of homogeneous SL(2, IR) metrics. This class does not admit Einstein metrics. This 
basin goes degenerate, yielding for a pancake degeneracy whereby a two-dimensional geometry 
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survives: two of the components of the metric increase without bound while the other shrinks 
to zero. 

(iv) : The Heisenberg-Basin. It contains all three- metrics flowing towards the three-parameter 
family of homogeneous Nil-metrics. Again, this class does not contain any Einstein metrics. 
This basin too undergoes a pancake degeneracy. 

(v) : The 22(1, 1) -Basin, where 22(1, 1) is the group of isometries of the plane with flat Lorentz 
metric. It contains all three- metrics flowing towards the three-parameter family of homogeneous 
Solv-metrics. Also this basin fails to contain Einstein metrics. This basin eventually exhibits 
a cigar degeneracy: the curvature dies away, and while one diameter expands without bound, 
the other two diameters shrink to zero. 

(vi) : The E(2)-Basin, where 22(2) is the group of isometries of the euclidean plane. It 
contains all three-metrics flowing towards the three-parameters family of homogeneous Solv- 
metrics containing the flat geometry. This basin is eventually attracted by flat metrics. 

(vii) : The 22 (3) -Basin, where 22(3) is the group of isometries of hyperbolic three-space. It 
contains all three-metrics flowing towards the one-parameter family of homogeneous metrics 
constant multiples of the standard hyperbolic metric. This basin is attracted to hyperbolic 
space. 

(viii) : The SO (3) x IR 1 -Basin. It contains all three-metrics flowing towards the two- 
parameter family of homogeneous metrics obtained by rescaling the standard product metric 
on S 2 x IR 1 . It does not contain Einstein metrics. This is a singular basin, it is attracted 
towards a curvature singularity: the round two-sphere shrinks while the scale on IR 1 , (or if you 
prefer, the S 1 factor in the original manifold), expands. 

(ix) : The H 2 x IR 1 -Basin, where 22(2) is the group of isometries of the hyperbolic plane. 
It contains all three-metrics flowing towards the two-parameter family of homogeneous metrics 
obtained by rescaling the product metric on the product manifold, IR 1 x H 2 . Again, this basin 
does not contain Einstein manifolds, and it is attracted towards a pancake degeneracy. 

The basins of attraction just described, and in particular those yielding for fixed points 
(Einstein manifolds), are relatively uninteresting in connection with the Renormalization Group 
interpretation of (|5j). Such fixed points, e.g. the round three-sphere, or the flat three-tori, 
are all totally attractive. As already recalled they can be thought of as distinct stable phases 
yielding for distinct cosmological models, characterized by the nature of the metric (and of 
its infinitesimal deformations) at the fixed point. For instance, the S'L r (2)-Basin characterizes 
FLRW models with closed spatial sections, and the related homogeneous anisotropic models, 
(see section 6 for more details). 
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4.3 Critical fixed points: an example 

The existence of critical fixed points, (critical in the sense of (]62])), characterizing the (universal) 
properties of (continuous) phase transitions between two different cosmological regimes, cannot 
be immediately read off from the above analysis. It is rather the consequence of the (conjec- 
tured) existence of the decomposition of a manifold into pieces that is associated, according to 
Hamilton, to the local singularities of the flow (^), (see the second step of Hamilton-Thurston 
geometrization program). The origin of this connection between critical behavior of (|52|), in 
the sense of Renormalization Group, and Hamilton's program can be seen by considering the 
following detailed example. 

Let us assume that topologically E is a three-sphere, E ~ S 3 . We shall consider on E a 
metric g\ obtained by glueing through a smooth connected sum two copies of a round three- 
sphere 

E-Sf^), (63) 

and endowing each S%y i = 1,2, factor with a round metric of volume v\ = 1, and the joining 
tube 

S 2 x ([0,1] C 1R 1 ), (64) 

with the standard product metric of volume const. e _3T , (r is a suitable parameter, see below), 
for r > 1. 



To explicitly construct this latter metric we can proceed as follows: 

Let yi and y 2 respectively denote two chosen points in both factor copies, and S^ 2 y 
Let h^y. 1R 3 — * Sf^, i = 1,2, be two imbeddings given by the exponential mappings 

exp yi : T^-ffi. 3 -^, 

exp y2 : T y2 Sf 2) ^JR 3 ^Sf 2) . (65) 
We assume that hr\\ preserves the orientation while h(2) reverses it. 

Let a: (0, oo) — > (0, oo) denote an orientation reversing diffeomorphism, and define a^: ]R 3 /{0} 
IR 3 /{0} by the map as(v) = a(\v\)-^ for every vector v G K 3 /{0}. For every point x\ = h\{v) 

in the geodesic ball B(y\,r)/{0} C 5 3 ^, with < r < ^, (with tt/2 the injectivity radius of 
the unit three-sphere S^), we identify hi(v) with h-zip-siy)) G B(y 2 ,r)/{0} C S^ 2 y 

The space obtained in this way 

Sf^Sf 2) = (Sf 1} - { yi }) U h(2)Q3h - 5 (5 3 2) - {y 2 }) (66) 

is a particular realization of the connected sum J30| of two copies of unit three-spheres. 



In order to give to the neck, joining the two Sf^, a cylindrical shape we blow up [17] 
the metrics of the three-spheres in the neighborhoods of the points y\ and y%. Consider, for 
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simplicity, only the factor since the argument goes in an analogous way also for the 
factor. 

Exploiting the exponential mapping the round metric of S 1 ^, (actually any sufficiently 
smooth metric), can be written in a neighborhood of y\ in geodesic polar coordinates as 

g{x) = dr 2 + r 2 h ij de i d9 j + 0{r 4 ), (67) 

where r = dist(x,yi) is the distance between y\ and the point considered, hijd9 l d9 3 is the 
metric on the two-dimensional unit sphere S 2 , and as usual, the higher order correction terms 
involve the curvature. 

Now, if we blow up this metric by rescaling it through r 2 we get, up to curvature corrections, 
a distorted cylindrical metric 

g( x ) = i^ = ^- + hlj de i d6 j + 0(r 2 ). (68) 



In order to eliminate the axial distortion due to we substitute for the radial variable r, 
introducing a new coordinate r defined as 

r = exp[— t]. (69) 

When expressed in terms of r we get that the blown up metric g(x) reduces to 

g(x) = dr 2 + h ij d6 i d9 j + 0(e- 2r ). (70) 

Thus, the blown up metric approaches the cylindrical metric exponentially fast as r — > oo. 
In order to introduce smoothly such a metric on the neck of S^ftS^ we can proceed as follows. 

Choose smooth functions 5^(r), i = 1,2, satisfying 

1 7T 

<%)( r *) = 12 > 0<r* < -, 

i 

= C°° , decreasing, ^ < ri < | , (71) 

7T 

<5(j)(r;) = 1, T-i > -, 

where = dist(x, y^)- Also, let us introduce the following inclusion maps 

X(D = ^i)-{w}^(^i)-{lft})U fc(!oa , fc -i(Sf 2) -{ltt}), 

X(2) : S? 2) - {y 2 } - (Sf 1} - { yi }) U fc(a)asfc -i (Sf 2) - {y 2 }), (72) 

X( 3 ) = M^Sfi) - {0}) U h 2 (T y2 Sf 2) - {0}) ^ (5^ - {yi}) U h(2)Qsh -i (Sf 2) - {y 2 }), 



26 



{X(i) and X(2) denote the inclusions of the spheres, minus the points yuy into the connected 
sum; X(3) is the inclusion of the neck). 

If {ia} is a partition of unity associated with the covering corresponding to the above 
inclusions X(a)> with & = 1)2,3, we can then define the following metric on S^sftS^y 

9{x) =^2taXl a )[g{x) ■ %)(dist(x,y (i) ))]. (73) 

a 



For x not in the geodesic balls (of radius vr/2), centered on y\ and ?/2, this is the standard 
round metric on the unit three-sphere; for x in the geodesic balls of radius 7r/4, centered on 
yi and 2/2 > this is, up to curvature corrections, the cylindrical metric introduced above. For x 
in the annuli between such balls g is a smooth interpolating metric joining the spheres to the 
cylindrical neck. 

By expressing ( [73]) in terms of the variable r, we get the metric on SqAS? 2 \ which is the 
round metric on each S^/B(yi, e~ T ) factor, and these factors are connected, for r large enough 
(i.e., nearby y^), by a thin flat cylinder. 



The Ricci-Hamilton evolution of (S%^S9 z \,g) can be explicitly constructed as follows. 
According to [25] let us write the metric on S 2 x 1R 1 as 



fflneckof Sf^Sf^ = D 9jr 1 + E 9s 2 > ( 74 ) 

where g^i is the metric on IR , g S 2 is the round metric on the unit two-sphere, and D and E 
are constants. The Ricci-Hamilton flow (|5^) preserves the structure of this metric and reduces 
to the coupled system of ordinary differential equations 

-E = --, 

drj 3' 

i° - m 

which immediately integrate to 



drj 3\E 



E = E -~r), 

D = ^> (76) 

where E 2 = E 2 (r/ = 0) is the initial radius of the round two-sphere, whereas Dq = D(r] = 0) is 
the scale on IR 1 . 
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Given this solution of (52), we can construct the Ricci-Hamilton evolution of (S^AS^yd) 
by looking for a solution of (|52"|) in the form 



(77) 



g(x;i]) = Y,taX* a )(v)[g(x) ' %)(dist 7? (x,y (i) ))], 



where the inclusion maps \{a) depend now on the deformation parameter r]. We assume that 
the metric g(x), (the round metric on the spheres and the standard product metric on the 
neck~ S 2 x 1R 1 ), being locally homogeneous, is preserved by (52) since the Ricci-Hamilton flow 
preserves isometries. As the plumbing between the spheres and the neck shrinks as r] increases, 
(as the S 2 factor in the neck), the inclusions X(a) are necessarily ^-dependent. The dynamics 
of X(a) can De obtained as follows. 

The Ricci-Hamilton flow ([52]) for g(x; rj) is given by 



dg a b(v) 

drj ~ 3 

9ab(v = 0) = g ab , 



2 < R( V ) > g ab (rj) - 2R ab (rj) 



and, in terms of the pulled-back metric it is 

d_ 

drj 



- 2R(x* {a) ^)g(x) ■ 5 {l) ) 



ab ■ 



(78) 



The left hand side of (|78[) can be evaluated according to a suggestion first exploited by D. 
DeTurck, viz., 



d r d 



drj 



+ X* a )(v) L w( V )[9 • 5(i)]ab(x{V,x)) 



(79) 



where the Lie derivative L w ^ [g ■ 5^} ab (x{rj > x )) is along the vector field w(n; a) which generates 
the ^-evolution of the inclusions X(a): ™ z -i 



dX(a)(v) 

drj 



w{V,X(a){v)), 



X(a)(V = fy = X(a)- 



(80) 



Let us denote by < R >( a ) the average of the scalar curvature R(x( a ){ x )) over the images 
of the inclusions X(a)> (viz., for a = 1,2, < R >( a ) is the average over S%y while for a = 3, 
the average is over the neck). In terms of these averages, the Ricci-Hamilton flow ( |78[) can be 
written as 
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d_ 

drj 



[<%) • 9ab(X(a){x))} 



= 3 < R(V) >(a) <%) • ff a6 (X(a)(^)) - 2 J R a fe(X(a)(^)) 

2 r 

+ 3%) • 9 ab [X(a) W) [< > - < #W >(o) 

- itu(r 7 )b-%)]ob(x( r ?) :c ))- 



(81) 



As can be checked, the Ricci-Hamilton flow preserves the local homogeneous geometries 



over the spheres S^, and on the neck, if and only if the vector field 10(77; a ) satisfies 
2 r 1 

0%) ■ 9 ab (X(a){x)) < R(rj) > - < R(rj) > {a ) = L w{v) [g • 6 {i) ] ab {x{v,x)). 



(82) 



By taking the trace of (82) with respect to g ■ Su\, we get 



< R> - < R >( a )= V k w k (r/; a) 



(83) 



where V denotes the riemannian connection with respect to g ■ 5 



■°(0- 



For a = 1,2, i.e. for the punctured spheres — {yi}, the above relation yields upon 



integration over — B(y, 



<R>- <R> 



(a) 



vol (SL - B(y h r)) = / w k (r); a)da k 



(84) 



where S (rf) is the //-dependent boundary of — B(yi,r). Since < R > — < R >(«), for 
a = 1,2, is proportional to the average scalar curvature of the neck, the term on the left 
hand side of (84) blows up as {Eq — |??) _1 as 77 increases. On the other hand, by the S 2 {rf)- 
rotational symmetry, the term Wk(f];a)da k is spatially constant on the two-sphere boundary 
S 2 (rj) of the punctured three-spheres. Thus, flS4"D implies that as r\ increases, the surface 
area of S 2 {r\) shrinks to zero, as {Eq — |r/) _1 , by moving along the outer normal direction of 
S 2 { V )c{Sf i) -B{y i ,r{r,))). 

For a = 3, i.e., for the neck, fl84|) yields 



<R> - <R> 



(a=3) 



vol (S 2 {rj) x IR 1 



w k {ri;a)da k - I w k {rj;a)da k , (85) 



where S 2 2 ^{rj) — S 2 ^{r]) is the oriented boundary of the cylinder S 2 {rj) x IR . Since < R > — < 



R >( a =3) is, up to small correction terms (coming from the collars joining the spheres to the 
neck), the average curvature over the spheres 5?a, % = 1,2, ( |85| ) simply tells us that as the 

S 2 {rj) factor in the neck shrinks, the scale-length of the IR 1 factor grows, so that the volume 
vol (S 2 ^) x IR 1 ) remains constant during the Ricci-Hamilton evolution, namely, 
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J S 2 i) w k (i 1 ;a)da k -J S 2 {) w k (r } ;a)da k 
vol(S 2 (r])xM^ - (2) (1) 



< R > - < R >( a=3 ) 



Notice that by introducing, as above, a new variable 

dist v (x,y(i)) = exp[-r(r/)], (87) 

with 

<n) = tt-Vt' ( 88 ) 

and by rescaling the unit two-sphere metric h{j according to 

hiM = [ l -\'n\hij, (89) 



the above analysis of the Ricci-Hamilton flow ( |7q ) provides on the neck the metric 

g(x; rj) = dT 2 ( V ) + h ij {r ] )d9 i d9 j + ©(e" 2 ^), (90) 
and on Sh, the standard round metric. 

Notice also that ( |90|) goes cylindrical exponentially fast in r\. As r\ increases < r] < |, the 
neck becomes longer and longer while getting thinner and thinner. In the limit, r\ — > |, we get 

iSysf 2) ,~g)^Sf 1) ]JSf 2) , (91) 

where the three-spheres carry the round metric of volume 1, while the smooth joining 
regions shrink exponentially fast around the points y^, i = 1, 2. 



4.4 Critical surfaces and topological crossover 

Let g Q be a metric on £ with positive Ricci curvature, 

Ric(g ) > (92) 

and with volume vol(T,,g ) = 2. 

By means of g Q and the metric g introduced in the previous section, we can construct on 
S a smooth one-parameter, (0 < t < 1), family of metrics gt, with g t= o = g and gt=\ = g by 
setting 

g t = (1 - t)g + tg, 

0<t<l. (93) 
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According to Hamilton's theorem p2}| , there is a right-open neighborhood of t = such 
that all metrics <ft in this neighborhood are in S'C/(2)-Basins and are attracted, under the 
action of the Ricci-Hamilton flow (|52|), towards the round metrics on 5 3 with volumes vt, 
(since the volume of gt changes with t, we have a family of fixed points, namely, round three- 
spheres parameterized by the corresponding volumes v% which are kept constant under the 
Ricci-Hamilton flow; thus Vt=o = 2). 

On the other hand, according to the remarks above, we have an open neighborhood of 
t = 1 such that all metrics gt in this neighborhood go singular under the Ricci-Hamilton flow. 
Indeed, for t = 1, the Ricci-Hamilton flow of g fixes the round S 3 factors while the joining 
tube, S 2 x IR , is driven towards a curvature singularity. By continuity, this behavior extends 
to a left-open neighborhood of t = 1, and a three-sphere in this neighborhood splits apart into 
two round spheres. 

It follows that there is a neighborhood of the g metric such that some of the metrics in this 
neighborhood are driven towards attractive SU (2)-Basins, while others are driven towards a 
singular (5^)J5^)-Basin. 

The set of metrics driven towards this basin defines a critical surface. Below we shall define 
a critical fixed point. Since the Ricci-Hamilton flow preserves the metric g, up to a trivial 
rescaling of the neck in S?i\$S%\, (rigorously speaking this is true only up to exponentially 
small correction terms which arise in the regions joining the three-spheres with the S 2 x IR 1 
neck), we can characterize the critical fixed point for (|52|), in the example considered, as the 
metric on •S^tt^) Si yen by 

3 

Gcrit = {g(x;v)\ V -> Tj}' ( 94 ) 

where g(x;n) is the one parameter //-family of metric solution of the Ricci-Hamilton initial 
value problem (^). This is a consistent characterization of a critical fixed point for (|52[), 
since (|52]) fixes (S^AS^, Gcrit), and as required for (]62]), the correlation length associated 
with (S^ftS^, Gcrit) is L = oo. This is so since, roughly speaking, in order to describe 
(5 , ^i)tJ5'^2)' Gcrit) we fix any representative (S^USj^) which is characterized by a correlation 
length L ~ the length of the neck, and then rescale with the Ricci-Hamilton flow (|78|) from L Q 
to L = oo. 



Notice however that in this particular case, the critical point is not to be related to a 
phase transition. As mentioned before, we are here in presence of significant finite-size effects 
which are concerned with a dimensional crossover and they show up, as usual, as an effective 
reduction of dimensionality In this case, the three-dimensional neck of (Sfx)$S?2\) goes 
one-dimensional. This crossover to (quasi-)one-dimensional behavior is not accompanied by a 
singular behavior in thermodynamical quantities, such as correlation functions, but nonetheless 
an anomalous behavior is present. Indeed, geometrically speaking, the critical fixed point 
(S^US^, Gcrit), and the corresponding critical surface, separates two stable phases under the 
renormalization generated by (52). One is given by the manifolds eventually evolving towards 
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the round three-sphere of volume 2. The other is generated by those manifolds which pinch 
off and eventually evolve towards two round three-spheres, each of volume 1. The pinching off 
through thinner and thinner necking is necessary for such a topological crossover. 



From a physical point of view, and as we shall see in section 6, the fixed point described 
above separates two possible different closed FLRW regimes. One with closed spatial sections 
which, at some particular instant, are a three-sphere of volume 2, while in the other regime we 
have two distinct closed FLRW universes having spatial sections (at a given instant) of volume 
1. We may have also many different necks corresponding to a regime whereby the spatial 
section £ yields for many closed FLRW universes. 



It is clear that the above explicit construction of a critical fixed point for (52) can in 
principle be generalized to more general situations. The strategy is to take two or more trivial 
fixed points, such as those associated with the SU(2) and -ff (3)-Basins, and connect them 
through the degenerate basins (such as SO (3) x IR 1 , as in the above example, or through the 
H(2) x 1R 1 -Basin, etc.). In this connection notice that the connected sum mechanism yielding 
for the (S^ftS^^-critical fixed point, can be generalized to an operation of joining two (or 
more) manifolds (corresponding to stable attractors) along tubular neighborhoods of surfaces 
(rather than points, as in the case for the standard connected sum). 



A particularly interesting connecting geometry would be H (2) x IR , (by compactifying the 
hyperbolic plane in a closed surface). In this latter case the scale of the hyperbolic geometry 
goes to infinity under the Ricci-Hamilton flow (pancake degeneracy p5|). Finite size effects 
are again at work, but this time the effective reduction of dimensionality is more interesting 
than the one in the (5'( 3 1 )tlS'( 3 2)> Gent) case. Indeed, two out of three dimensions are infinite, 
and there will be a crossover to a critical behavior with critical exponents characteristic of a 
two-dimensional system. 



One can consider the analysis presented above as a physically non-trivial application of 
Hamilton-Thurston's geometrization program. In some of its aspects it is rather conjectural 
and speculative, but in our opinion, it is quite intriguing that motivations coming from geometry 
and a physical problem, like the one addressed here, namely, the construction of cosmological 
models out of a local gravitational theory, go hand in hand in such a way. 



5 Linearized RG flow 



The relative slopes of < / > tm (g±) and < / > tm (#2), with / given by equation (35), as 
m —* 00, and for g\ in a neighborhood of g<i > are of some relevance to our discussion. In a 
standard RG analysis such relative slopes are related to critical exponents. Given the blocking 
procedure for < / > €m (g\) as m — > 00, yielding for a g± "renormalized" according to the Ricci- 
Hamilton flow (52), one can sensibly ask what happens if g\ is slightly perturbed, namely, if 
we replace it by 

9ab -> 9ab + 5K ab , (95) 
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where K ab is a symmetric bilinear form (a choice of the symbol is quite intentional, since later 
the above consideration will be applied to the second fundamental form). It can be shown that 
if g\ is scaled, according to the Ricci-Hamilton flow (|52|), then K ab gets renormalized according 
to the linearized Ricci-Hamilton flow, namely (rj in the brackets suppressed), 

^-K ab = l < R > K ab + lg ab {\<Rg ab K ab >-\<R><g ab K ab >- 
or\ 3 3 2 2 

< R ab K ab >] - A L K ab + 2[div*(div(K - ^(Tr K) g))] ab , (96) 

with the initial data K ab (r\ = 0) = K ab , where K 6 S 2 T, is a given symmetric bilinear form, 
Ai is the Lichnerowicz-DeRham Laplacian on bilinear forms 

& L K ab = -V s V s K ab + R as K% + R bs K s a - 2R asbt K s \ (97) 

and the operators A^, div* , div and TV are considered with respect to the flow of metric 
(fiS 7 ?) ~~ * 9{f])i solution of (|52[). The div (here, minus the usual divergence) is the divergence 
operator on 5 2 S, div* is the L 2 - adjoint of div , acting from the space of vector fields on E to 
S 2 T, (it can be identified with ^ [Lie derivative] of the metric tensor along a vector field). 

Note that K(rj) solution of the linear (weakly) parabolic initial value problem ( |96|) always 
exists and is unique |22[| , and represents an infinitesimal deformation of metrics connecting the 
two neighbouring flows of metrics, g(rj) and </(??), (obtained as solutions of problem (|52|) with 
initial data g(rj = 0) = g and g'(rj = 0) = g(rj = 0) + eK(rj = 0) + C(e 2 ), respectively). For 
what concerns the structure of this solution, one can verify that corresponding to the "trivial" 
initial datum K(r\ = 0) = Lxg (where X : E — > TE is a smooth vector field on E), the solution 
of (^) is provided by 

K a M = L x g a b(v)- (98) 

This property expresses the Di//(E) equivariance of the Ricci-Hamilton flow. (Notice that X 
is ^-independent). 

The above fact follows by noticing that along the trajectories of the flow (r],g) — > g(rj), 
solution of (|52[), we have 

8 2 

-Q^L x gab{v) = Lx[-r^9ab(r])} = g < R(v) >v L xg a b{v) ~ ^L x R a b(v)- (99) 



But the Di//(E)-equivariance of the Ricci tensor, i.e. the fact that Ric((p*g) = cp*Ric(g) for 
any smooth diffeomorphism ip : E — * E, implies that 

L x Rab = D Ric(g) ■ L x gab, (100) 
where D Ric{g)K is the formal linearization of Ric(g), around g, in the direction K: 

D Ric(g) ■ K = -[Ric(g + tK)} t =o 

= -A L K - div*[div(K - -(TrK)g)]. (101) 
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Upon introducing ( |100| ) in (99) we get 
2 

Q^Lxg a b(ri) = -< R(v) > v L x9ab(v) ~ 2D Ric(g(rj)) ■ L x g a b(v)- (102) 

One can check that the right hand side of the above expression coincides with the right hand side 
of (96), when this latter is evaluated for K ab (r]) = Lxgabiv)- Hence Lxgabi 7 ]) solves the partial 
differential equation ( p6| ) and, since for rj = 0, K ab = Lxg a b, the uniqueness of any solution of 
the initial value problem (|96|) implies that K ab (rj) = Lxgab(v)> whenever K ab (r) = 0) = Lxg a b, 
as stated. 

Moreover, if K (rj) is a solution of (|96|), with initial datum K(r] = 0) = K, then the space 
average of Tr K(rj) over (T,,g(rj)) is preserved along the flow (r],g) — * g(rj), namely, 

< Tr K(rj) > v =< TrK > a , < r? < oo. (103) 



This property of the solutions of ( |9q ) is an immediate consequence of the volume-preserving 
character of the Ricci-Hamilton flow. 

Finally, another relevant property of the initial value problem (|9^) can be stated as follows. 
If (rj,K a b) — > K a \,(ri) is the flow solution of (|96|) , with initial datum K a b(r] = 0) = K a b, then it 
can always be written as |35| 

K a b{ri) = K ab (ri) + L v ( v) g ab (r]), (104) 

where the bilinear form K ab (rj) and the 77— dependent vector field v(rj), respectively, are the 
solutions of the initial value problems: 

^-K ab = l < R > K ab + \g ab [\ < Rg ab K ab > -\ < R >< g ab K ab > - 
or/ 3 3 2 2 

< R ab K ab >] - A L K ab , 

ka^ri = 0) = K ab , (105) 



and 



d 1 

Q-v a {ri) = -V c {K ca - -k rs g rs g ca ), v(rj = 0) = 0. (106) 



To summarize, as rj — > 00, K ab (rf) may either approach a Lie derivative term, such as 
L v ^g ab (r]), or a non-trivial deformation k ab (rj) |^5| . The non-trivial deformation is present 
only if the corresponding Ricci-Hamilton flow for g a b(v) approach an Einstein metric on X 
which is not isolated. In such a case, (e.g. flat tori), there is a finite dimensional set of such 
Einstein metrics, and the non-trivial k ab simply are the infinitesimal deformations connecting 
one Einstein metric g\ in X and the infinitesimally non-equivalent one. Also in this case the 
Lie derivative term may be present. What it represents is a reparameterization of the metric 
gi (under the action of the infinitesimal group of diffeomorphisms generated by v (see equation 
(|104j )), ("gauge artifact"). This latter Lie derivative term is the only surviving term when g is 
isolated (like e.g. in the case of the round three-sphere. As is known, the round metric g on 
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the three-sphere S s is isolated, in the sense that there are not volume-preserving infinitesimal 
deformations of g mapping it to another inequivalent constant curvature metric g 1 '.) 

To summarize, the flow fll05|) , ( |106| ), tells us how to renormalize the second fundamental 
form in such a way that the blocking prescription < / > €m — > < / > £m+1 works both for the 
initial metric g as well as for the perturbed metric g + bK ab . 



5.1 Scaling and critical exponents 



From the characterization of critical fixed points g* b for (52), (see section 4.2), we can get 
information on the critical exponent characterizing critical behavior of the metrics nearby g* b . 
This is the content of this section. In particular, we shall discuss the critical exponents related 
to the critical fixed point (SfufySSn, Gcrit)- Even if this point is not a thermodynamically 
interesting critical point, it exhibits many of the general features of the more interesting type 
of singularities. 

In order to characterize these critical exponents we can use the linearized Ricci-Hamilton 
flow associated with the one-parameter family of metrics (,S%\$S%\-, Gcrit)- I n the following we 
shall however proceed more directly and examine the properties of the two-point correlation 
function associated with {S%AS%y Gcrit) and the probability law of relevance to our analysis. 

Let yi £ Sf^, i = 1,2, be the two points in (S^^S^, Gcrit) around which the round metrics 
of S%-> have been blown up. Let / denote a non-negative scalar field on (S^ftS^), distributed 
according to the probability law formally defined by 

e W [-H(f)]U x df(x) 

Jexp[-H(f)]U x df(xY [ } 

where the functional integration is over the space of fields /: {Sfi^S^) — ► M, equipped with 
the L 2 -inner product 

(/I/O = I , f(x)f'(x)d fJ , g (x), (108) 

and where 

H(f) = f f{x)d N {x). (109) 

(We could consider f as related to the matter fields, as in / = ap + a 1 Ji, but the following 
analysis is quite independent from a particular meaning of /). Notice also that the measure 
dfi g (x) in the above formulae is the riemannian measure on (5'( 3 1 )tt'S'^) associated with the 
metric g{x) defined by (fT5|). 



Let us now concentrate on the behavior of ( 107 ) when the neck of (S^^S^) gets thinner 
and longer under the action of (|52|). 

It can then be checked that along the Ricci-Hamilton flow ([76]) associated with g| nec k, 
both the L 2 -inner product ( |108|) and the Gibbs factor exp[—H(f)] corresponding to ( 109 ), are 
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invariant. Thus it follows that the probability measure ( |107| ) is invariant under (52), and the 
correlation function defined by 



E d p[f(yi)f(y 2 )}, (110) 

where E d p denotes the expectation with respect to dP, is well defined over the critical fixed 
point Q crit . 



In section 4.1 we have interpreted the Ricci-Hamilton flow ( p2| ) as the RG flow in a finite 
geometry, characterized by the length scale L, which is large compared to the microscopic scale 
(in particular, it is much larger than the radius of the geodesic ball coverings used to discretize 
the theory). The correlation function depends on such a dimensional parameter. If 

L( yi ,y 2 ) = / VDdz (111) 
Jyi 

denotes the distance between the points y\ and y 2 along the cylindrical neck of (<S'( 3 1 )tt5^2))' 
and if this distance is large (as compared with the radius of geodesic ball coverings), then the 
correlation length £ associated with the two-point (connected) correlation function can be read 
off from 

£~ ~ L (^'^) ni „x 

* \uE dP {f{y 1 )f{y 2 )) conn 

Since the correlation function remains invariant under the Ricci-Hamilton deformation of the 
cylindrical neck, we get that along (|76j), £/L(yi,y 2 ) remains constant which implies that on 
0Sm\ tlt?(2) > Qcrit) the correlation length £ behaves as 

£~ y/D^Eo[Eo-lrj\- W . (113) 



According to standard usage, we can define the critical exponent v associated with the 
correlation length of a finite size system (with typical size L) by the condition 

^(L,P)\ P = PC ^L 1+ K (114) 

where p is a parameter driving the system to criticality, and p c is its corresponding critical 
value. In our case, it is natural to set p = rj, with p c = rj c = 1 . This immediately yields for the 
critical exponent u the value 

v = l. (115) 



A similar computation for the critical exponent associated with the correlation length can 
be carried out when the connecting geometry is H 2 x IR 1 . This takes place when two riemannian 
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manifolds Mi and M2, which are supposed to evolve nicely under the Ricci-Hamilton flow, are 
connected through a tubular neighborhood of a surface Sh of genus h, viz., a = M\$s h Mi. 



In this case, the metric on the neck can be written as [25| 



Qneck = Dg u i + Eg H 2, (116) 

where 17^-2 is the metric on the hyperbolic plane. The Ricci-Hamilton flow equations take a 
form similar to the S 2 x IR 1 case, up to a minus (important!) sign, namely, 

d -E 2 



drj 3 

—D = - 
drj 3 \E 



d 4 (D\ 

D = -o l-B > ( 117 ) 



which upon integration provide 



2 

E = E + -7], 

where E 2 = E 2 (j] = 0) is the initial scale of the hyperbolic geometry, whereas -Do = D(tj = 0) 
is the scale on ]R . 

The scale E of the hyperbolic geometry increases linearly with 77, while the scale factor D 
on the line IR 1 shrinks. It can be checked p5j ], that corresponding to this scale dynamics, the 
curvature decays according to ||-Ric|| = \^2/(Eq + |r/), and we get in the limit 77 — » 00 a pancake 
degeneracy. 

The relevant correlation function is now Edp(f (yi) f (y2))conn, with y\ and j/2 fixed points 
in the H 2 factor, (i.e., on the surface Sh). Again, owing to the symmetries of the geometry 
involved, it immediately follows that the ratio between the correlation length £ and the distance 
£(2/1,2/2) must remain constant under the Ricci-Hamilton flow. This implies that the correlation 
length behaves as 



, /••() • I' I- (H9) 



to which corresponds again a critical exponent v = 1, (in this case one cannot apply ( |114 ) since 
the system is actually going to an infinite size). 

A striking feature of these topological crossover phenomena, associated with the renor- 
malization of the cosmological matter distribution, is that their pattern resembles the linear 
sheet-like (or sponge- like) structure in the distribution of galaxies on large scales (e.g. [pi). It 
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is evident that if, (E,g, K, p, J), the initial data set for the real universe (see the next section), 
is close to criticality, in the sense discussed in the previous section, then the corresponding 
averaged model exhibits a tendency to topological crossover in various regions (the ones where 
the inhomogeneities are larger). Filament-like and sheet-like structures would emerge, and 
the overall situation would be the one where such structures appear together with regions of 
high homogeneity and isotropy, in some sort of hierarchy. This situation is akin to that of a 
ferromagnet nearby its critical temperature, whereby we have islands of spins up and down in 
some sort of nested pattern. Even if there is a tendency to homogeneity at very large scales, 
the picture just sketched, of the "hierarchy of structures", might be qualitatively valid in a 
good part of the universe; indeed recent observational data seem to suggest the existence of 
still larger and larger structures (e.g. [13|). Notice that this picture bears some resemblance to 
the "cascade of fluctuations" in critical phenomena |^9| . Droplet fluctuations nucleated at the 
lattice scale in the critical state can grow to the size of the correlation length where the details 
of the lattice structure become lost and the scale invariant distributions of the large "droplets" 
are universal. 

A possible (though not yet clarified) connection of this whole picture with the self-organized 
criticality (SOC) |3j can be envisaged, whereby the problem of structure formation, in terms 
of growth phenomena, could be tackled in the framework of avalanche activity used in SOC 
(cf. ^(|). It would be particularly interesting to estimate the scales in the universe, where it 
is necessarily critical and trapped into self-organized (critical) states (similar suggestion was 
recently posed in p^j ). 



6 Effective cosmological models 

The results of the previous sections have interesting consequences when applied in a cosmolog- 
ical setting. 

Let us assume that at the scale over which General Relativity is experimentally verified, 
a cosmological model of our universe is provided by evolving a set of consistent initial data 
(X, g, K, p, J), according to the evolutive part of Einstein's equations. 

The data (£, g, K, p, J) describing the interaction between the actual distribution of sources 
and the inhomogeneous geometry of the physical space, (X,g,K), are required to satisfy the 
Hamiltonian and the divergence constraints. Respectively, 

8ttG P = R(g) - K ab K ab + k 2 , (120) 

V i K lh -V h k = 16irGJ h , (121) 

where, k = K a a , (see also section 4, where they were written in terms of the three-metric g a b 
and its associated conjugate momentum 7r a b). 

According to the blocking and renormalization procedure, discussed at length in this paper, 
we can implement a coarse-graining transformation on this actual data set by suitably renor- 
malizing the metric and the second fundamental form, while retaining the functional form and 
the validity of the constraints. 
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The metric is renormalized according to the Ricci-Hamilton flow (52), 

f %M =l<R(r ] )>g ab (r ] )-2R ab (rj) 
\ gabiv = 0) = gab- 



In order to renormalize the second fundamental form, let us rewrite it in terms of the 
deformation tensor H a b defined by 

H a b = K a b + L a g a b, (122) 

where a is the shift vector field providing the three- velocity of chosen instantaneous observers on 
£, with respect to which we are implementing the renormalization procedure. We renormalize 
K a b by rescaling this deformation tensor H a b according to (|9"6|), viz., 

^<R>H ab + ^g a b[\ < Rg ab H ab >-\<R>< g ab H ab > - 
< R ab H ab >} - A L H ab + 2[div*(div(H - ~(Tr H) g))U, (123) 



d_ 

drj 



H, 



ab 



with the initial data H a b(rj = 0) = (K a b + L a g a b)r]=o- Notice that we renormalize the defor- 
mation tensor H a b, rather than K a b directly, because in this way we can get rid of the possible 
L>i//-induced shear which may develop in lim^oo K a b(rj). Indeed, according to equation ( [104|) . 
as 7] — > oo the solution of this initial value problem, K a b(rj), approaches a non-trivial deforma- 
tion lim^oo K a b(rj) plus a Lie-derivative term, lim^oo L v ^g a b(rj). The former is present only 
if the Ricci-deformed metric g a b = hm^oo 9ab (v) 1S n °t isolated. 

Since at this stage we are mainly interested in FLRW space-times, let us assume that g a b 
is isolated, while in order to take care of the £>i//-induced shear, lim^oo L v ^g a b{ri), we can 
choose the shift vector field ot in such a way that lim^oo L v ^g a b(r]) is compensated by L a g a b- 
Since L a g a b is a trivial datum for (|96|), it is thus sufficient to choose 

a 1 = lim v(n). (124) 

Notice that ot is the three-velocity vector of the chosen instantaneous observers on S, 
thus ( |124| ) provides a map identifying corresponding points between the initial actual manifold 
(£,<7, K, p, J) and its renormalized counterpart hm r) ^ 00 (S,g, K, p, J)(rj). 

In this section, and mainly for actual computational purposes, we assume that the original 
inhomogeneous initial data set is such that the Ricci-Hamilton flow is global. As already 
mentioned, we are interested in connecting an inhomogeneous cosmological space-time to its 
corresponding FLRW model. This is the case, in particular, if we assume that the original 
manifold (£,<?) has a positive Ricci tensor, (this case is obviously quite similar to the analysis 
in |]], there are however important differences that we are going to emphasize). Or more 
generally, if we assume that the original manifold (T,,g) is in the SU ^2)-Basin of attraction or 
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nearby the critical point tti-S^, with i = 1, 2, . . ., yielding for a manifold (E, g), nucleating under 
the Ricci-Hamilton flow (|78|) , (extended to many connected sums), to disjoint three-spheres S^y 

Given this setting, we see that due to the properties of the Ricci-Hamilton flow we have 
lim^oo Kab(r)) = ■ g a b- The given K a b is deformed by gradual elimination of its shear 
K a b — \kgab and the original (position dependent) rate of volume expansion k is being replaced 
with its corresponding average value. 



Since the constraints, (120) and ( |121| ), are required to hold at each step of the renormal- 
ization procedure, we get 

BnG(ri)p(ri) = R(g(rj)) - K ab ( v )K ab ( v ) + k 2 ( V ), (125) 

V(r)yK ih (n) - V( V ) h k = ltorGfr) J h ( V ), (126) 
where we have explicitly introduced a possible 77-dependence into the gravitational coupling G. 



Let us now explore the consequences of ( |125| ) and ( |l26j ). From the stated hypothesis on the 
Ricci-Hamilton flow it follows immediately that, K a b(r]) — ► \kg a b as 77 —> 00, thus ( |126| ) implies 
that 

HmJ /l (7 ? ) = 0. (127) 



In order to analyze ( [1250 , we will make use of a property of the Ricci-Hamilton flow, namely 
that the flow K(rj), solution of (|96|), is such that -jL < k{rf) >r;= 0, i.e. the space average of the 
trace of the second fundamental form remains constant during the deformation. This allows us 
to write 

< k > 2 = k 2 , (128) 

since in the limit, the volume expansion is simply a constant, and where < . . . > denotes the 
full average of the enclosed quantity with respect to the initial metric. 



Equation (|12g| ) provides the Hubble constant on the FLRW time slice associated with the 
smoothed data corresponding to (X, g, K, p, J). 

More explicitly, let us write the FLRW metric in the standard form (units c = 1) 

ds 2 = -dt 2 + S 2 (t)da 2 , (129) 

where da 2 is the metric of a three-space of constant curvature and it is time independent. As 
we are interested in the three-sphere case, the metric da 2 can be written as 

da 2 = d X 2 + sin 2 X {d6 2 + sin 2 9d(p 2 ). (130) 

Since the volume vol(E, g) of the original inhomogeneous manifold is preserved by the Ricci- 
Hamilton flow, we can relate the factor S 2 (t), providing the inverse (sectional) curvature of the 
FLRW slice t = t a , to vol{E,g) 
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Notice in particular that the scalar curvature, towards which R{g(rj)) evolves under the Ricci- 
Hamilton flow, is given by 
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H.li, *<*,))- J,- i^m"'. (132, 



Having said this, the equation ( |125| ) becomes in the limit, after extracting a trace free part 
of K ab , 

8nGp = R + -P, (133) 

since no residual shear survives, and where we have introduced the renormalized gravitational 
coupling 

G = lim G(rj). (134) 

r\ — »oo 

On the other hand (|125|) gives 



-k 2 ^) = %irGp{rj) - R(g( V )) + k ab ( V )K ab ( V ), (135) 
where the shear K ab = K ab — ^kg ab has been explicitly introduced. 

We can rewrite the last equation upon taking the average with respect to the initial metric, 

as 

- < k 2 > = 8itG < p > Q - < R > a + < K ab K ab > . (136) 



By exploiting (|128 ) and ( |136D we can get the Hubble constant of the FLRW model we want to 
use to describe the time evolution of the renormalized initial data set lim„_ +00 (S,5,if, p, J)(rj). 

6.1 The renormalized Hubble constant 

Taking into account ( |128| ), and noticing that the Hubble constant on the FLRW slice corre- 
sponding, via the Ricci-Hamilton flow, to (E,g,K) is 

2 



we get 



Hi = ^ < p >o ~ < R >o +^ < K ab K ab > ~(< k 2 > -<k > 2 ). (138) 



This is a rather tautological rewriting of (137), justified by the fact that in this form H 2 clearly 



shows a contribution from the spatial average of the shear. In particular, if originally (i.e., for 
rj = 0) k is spatially constant, Q138 ) reduces to 



Hl = ^—<p> -\<R> +\< H ab H ab > +l< L a g ab L a g cd g ac g bd > D , (139) 
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where we have explicitly used the decomposition K ab = — L^^, (see equation ( 122Q ), and 



2 

L a g a b = V a a fe + V b a a - -g afe V c a c (140) 
is the conformal Lie derivative of the metric g ab along the vector field a. 
The above expressions for H Q correspond, through ( |133j ), to the matter density distribution 

8irGp = 8vrG < p > Q +R- <R> + < K ab K a b >o ~\{< k 2 > Q - < k > 2 Q ). (141) 



We have to consider this expression for Gp as the renormalized effective sources entering 
into the Friedmann equation, if we want to describe the real locally inhomogeneous universe 
through a corresponding idealized FLRW model. As expected, the renormalized matter density 
shows contributions of geometric origin, either coming from the shear anisotropies or from the 
local fluctuations in curvature and volume expansion. 

At this stage, a few comments are in order concerning, in particular, the effective Hubble 
constant H Q ( |13§| ), First, we wish to emphasize that this is the theoretical expression for the 
Hubble constant if one wishes to model the real locally anisotropic and inhomogeneous universe 
with a corresponding FLRW one. Expression ( |13§| ) clearly shows that apart from the expected 
contributions to the expansion rate coming from matter and curvature, there are two further 
contributions: 

(i) a negative contribution coming from the local fluctuations in the expansion rate. Apparently, 
this is not a very significant term, since we may choose, as a convenient initial data set describing 
the real locally inhomogeneous universe, a data set supported on a three-manifold S of constant 
extrinsic time, (i.e., k spatially constant on £). Actually, the optimal choice of the data set with 
respect to which smooth-out the lumpiness in matter and geometry is a deep and interesting 



question on which we comment later on, in section |6\3 . 



(ii) a positive contribution coming from the shear term | < H ab H a b > Q , and from the observer 
dependent term i< L a g a bL a g c dg ac g bd > Q which depends on the three- velocity a of the observers 
with respect to which the renormalization is carried out. 

These latter terms are the most important non-standard contributions to H 2 , and have 
their origin both in the presence of the gravitational radiation 

~< HfH Lab > G , (142) 

where H± is the divergence- free part of the deformation tensor H a b, and in the anisotropies 
generated by the motion of matter 

1 -7ab 



-<ff|fff|| a6>o , (143) 
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where Hu is the longitudinal part of the deformation tensor, obtained as a solution to the 
equation 

V\K\\) ih - \v h k = l6TrGJ h + V l (L a g lh ). (144) 



Recall that the shift vector field a is connected to the current density J by the requirement 
that it is chosen in such a way as to eliminate the longitudinal shear in lim^oo K ab (w). Also, 
notice that 

< H ab H ab > a = < HfH ±ab > a + < H^H ]]ab > o . (145) 

In practice, the deformation energy < H ab H ab > Q yields a contribution to the Hubble con- 
stant H Q which can be roughly estimated by exploiting the anisotropy measurements in the 
cosmic microwave background (CMB) radiation, as long as the frame used in averaging (i.e., 
the lapse a and the shift a 1 ) is, on the average, comoving with the cosmological fluid. We also 
require that the original locally inhomogeneous manifold (T,,g,K) does not differ too much 



from a standard FLRW t = const, slice. In such case one can apply the analysis of [36| to con- 
clude that, at the present epoch, the ratio between the deformation shear H ab and expansion 
is of the order 




(146) 



where e = max(ei, £2, £3) denotes the upper limit of currently observed anisotropy in the 
CMB radiation temperature variation, and where e\, 62, £3, respectively, denote the dipole, 
quadrupole, and octopole temperature anisotropies. On choosing e ~ 10~ 4 , as indicated by 
the recent CMB radiation anisotropy measurements, one gets that the shear deformation is at 
most about 1CT 3 of the expansion [^q] . 

Thus, as long as we assume that the original locally inhomogeneous manifold (T,,g,K) 
does not differ too much from a standard FLRW t = const, slice, the contribution to H a from 
< H ab H ab > is certainly quite small, and we can reliably write 

Ho ^ < P >o ~ < R >o L a g ab L a g cd g ac g bd > a . (147) 

The last term involving the shift vector 5 is a velocity effect term which is by no means 
small, at least a priori. Thus the renormalized Hubble constant needed, to describe by a 
model-FLRW a locally inhomogeneous and anisotropic universe, is not provided by a naive 
average even if the actual universe, to be smoothed-out in FLRW-modelling, is not too far from 
homogeneity and isotropy. Clearly, we are expecting that in such a case the actual contribution 
of the observer-velocity term is rather small, but this is a question which is difficult to handle. 
Technically speaking, we would need estimates on the size of the Dif f '-induced shear generated 
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upon smoothing by the Ricci-Hamilton flow. This is an issue under current investigation. We 
wish also to stress that in this velocity term there is hidden a non-trivial scale dependence. 
Indeed, assuming that we wish to model the actual universe by a FLRW one only up to a certain 
scale, then the shift vector a. needed to cure the Di//-induced shear, (which is tantamount to 
saying the proper selection of a frame with respect to which we are carrying out the partial 
smoothing), depends on such a scale, and the larger the scale the larger the contribution. 

There is another intriguing possible explanation for a larger than expected Hubble constant. 
Indeed, if we take seriously the possibility that the real universe may be close to the critical 
phase, as argued in the previous sections, then the contribution from the shear is not just of a 
conceptual value. For example, the original data set (T,,g,K) may be near the critical surface 
associated with the feSf* critical point. In this case we may generate, upon smoothing, a whole 
family of disconnected t = const. FLRW slices, each one with its own Hubble constant H (i). 
Explicitly, let us assume that the initial data set (£, g, K) is attracted upon averaging towards 
a critical (iS^fjS^tt . . . tt<S r L0 ) this being the case if the original metric g exhibits regions of 
large inhomogeneities. Assuming, for simplicity, that the rate of volume expansion is spatially 
constant, we get that each S%, factor, i = 1, ... ,n, inherits a Hubble constant provided by 

Hl(i) = ^<P>o-l< R >o+l< K ab K ah > . (148) 

These Hubble constants can be quite dominated by the large anisotropies < K^Kab > Q of the 
original manifold. Indeed, the previous estimate on the smallness of the shear term is, strictly 
speaking, valid only in the observable domain defined by our past light-cone from last scattering 
to the present day. Thus, if the factors resulting from the critical behavior of (S, g, K) are 
not smaller than the spatial sections intercepted by the interior of this light-cone, we would 
not notice the contribution from the local shear, (since we would have been looking at a rather 
homogeneous and isotropic island) , the large contribution would come from the regions of large 
inhomogeneities and anisotropies which, under the Ricci-Hamilton renormalization, undergo 
the topological crossover. Thus, a value of the Hubble constant may quite well depend on a 
possible large shear outside our observable domain. 

The behavior just described also shows that the RG procedure developed preserves the size 
of the regions where spatial homogeneity and isotropy is observed, even if large inhomogeneity 
may be globally present. This tells us what are the averaging scales on which a FLRW model 
is a good approximation. More precisely, in order to determine the averaging scale associated 
with a FLRW model we want to use, to describe (£,g,K), we need to understand more in 
detail the (conjectured) Hamilton-Thurston decomposition of (E, g, K). The specific example 
discussed in section |4J^ may serve as an indication. 

6.2 The scale dependence of the matter distribution 

In the above analysis, we also introduced a renormalized gravitational coupling. In a sense 
this is superfluous since the three-dimensional metric g of E is acting as the running coupling 
constant, and we can always reabsorb G into a definition of g. Nevertheless, the use of the 
renormalized coupling G may be helpful if one wishes to use the standard average of matter 
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< p > in the Friedmann equation, rather than the effective matter distribution p. The explicit 
expression for G can be easily obtained by setting p =< p > Q in (|141|) : 



8,0 - 8,0 + *- < R >° + < K K - k >° " >° dfl . (149) 

< P >o 

Notice however, that it is G{rj)p(rj) which is inferred from measurements for different scales, 
and thus the use of G is not particularly remarkable. 

In this connection, it is more important to discuss the dependence of G(rj)p(rj) as the local 
scale is varied, namely as r\ increases, (recall that rj is the logarithmic change of the cutoff 
length associated with the geodesic ball coverings). For simplicity, we do this only for the case 
in which no shear is present (K a b = 0), and the rate of volume expansion is spatially constant 
(k = const.). Under such hypothesis, we get for the scale-dependence of the average < Gp > 

^ <G ( V )p(r 1 )>] = |- < R(g( V )) >= 

= 2 < R ik R ik > +- (< R > 2 - < R 2 >) , (150) 



where R^ = Rik — ^QikR is the trace-free part of the Ricci tensor. From this expression 
we see that, not only shear anisotropies but, also metric anisotropies favour an increasing in 
G(n)p(r]) with the scale. To give an explicit example, let us consider as an initial metric to be 
smoothed-out a locally homogeneous and anisotropic S'C/(2)-metric g. Following the notation 
and the analysis in the paper of Isenberg and Jackson |2j| , we can write such a metric and its 
Ricci-Hamilton evolution, in terms of a left-invariant one-form basis {0 a }, a = 1, 2, 3, on SU (2) 
as 

g = AirjW 1 ) 2 + B(n)(8 2 ) 2 + C( V )(9 3 ) 2 , (151) 

where A, B, C are scale (r/)-dependent variables. With respect to this parameterization, the 
scalar curvature is given by 



R{rj) 



1 



[A 2 -(B- C) 2 ] + [B 2 -{A- C) 2 } + [C 2 -{A- B) 2 } 



While the squared trace-free part of the Ricci tensor is given by 

,2 I 

6 
1 



\\Ric\f 



[A 2 — {B — C) 2 } 2 + [B 2 -{A- C) 2 } 2 + [C 2 -{A- B) 2 } 2 
[A 2 - (B - C) 2 ][B 2 -(A- C) 2 ] + [A 2 — (B — C) 2 )[C 2 -(A-B) 



(152) 



(153) 



+ -[B 2 - (A - Cf][C 2 - (A 



B) 



The Ricci-Hamilton flow for this metric g exponentially converges to the fixed point A = 
B = C = 1, with the normalization ABC = 1, and with A(r]) > Birf) > C(rj) for all r\ [25]. 
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From the above expression for R(n) it follows that R(rj) monotonically increases from its initial 
value R(rj = 0) and exponentially approaches 

\imR(r,) = ^. (154) 

This increase is generated by the exponential damping of the anisotropic part of the Ricci tensor 



( 154 ) which is smoothed by the Ricci-Hamilton flow and, roughly speaking, is redistributed 
uniformly in the form of scalar curvature. 

Notice that if the anisotropy is large, the actual increase ^from R{n = 0) to lim^oo R(rj) = 
R may be quite significant. For instance, if for rj = 0, we have A(rj = 0) = 1, B{n = 0) = 2, 
C{n = 0) = 1/2, we get R{n = 0) = 7/8, while R = 3/2. The renormalized scalar curvature 
will have increased by nearly as much as 70% of its original value. 

Thus, if we smooth-out the initial data set (S ~ S s , g, K, p,J), with g the above SU{2)- 
metric, K a b = ^g a bk, with k = const., J = 0, and with the matter density p such that 
the Hamiltonian constraint holds, we get that G{rj)p{n) monotonically increases as n —* oo, 
exponentially approaching a fixed value. 

Under the above simplifying assumptions concerning the shear and the rate of volume 
expansion, we get 

8irGp = 8vrG < p > Q +R- < R > a (155) 

and, as remarked above, the correction term R— < R > G may be quite large, of the same order 
of magnitude as the naive average 8irG < p > Q . 

More generally, we can easily get an expression for the variation of G(rj) < p(rj) > as rj is 
increased from a given scale, say from i] to a larger scale rj = rj Q + 5n. 
From equation (|15C| ) we immediately get 

8vrG(r / ) < p( V ) > = 87rG(ry ) < p(rj ) > +2 < R k R lk > M 5n (156) 

+ - f < R > 2 - < R 2 >) 5r, + 0{5n 2 ). 
3 V ' (vo) 

With our simplifying assumptions, that there is no shear and that the rate of volume expansion 
is spatially constant, we can exploit the Hamiltonian constraint and rewrite (|157j) as 



G( V )<p( V )> = G{r ]o )<p{i lo )>+—<R ik R lk>{:no) 5n (157) 
- (rj ) < p(n ) > I — - >2 \5n + 0{5n). 

Thus, in leading order, a large density contrast <p ^^^y0^^ > a t a given observational 
scale r/ , tends to reduce the value of G(n) < p{rj) > as the scale of observation is increased. 
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This reduction effect should dominate on local scales where large inhomogeneities are present. 
On sufficiently large scales r) , cosmological data suggest that the density contrast tends to 
decrease, and in such cl CELS6, ELS the scale of observation is further increased, equation ( p.58|) 
shows that G{rj) < p{rj) > now tends to become larger, with deviations from the original value 
driven by anisotropies in the curvature. 

Thus, the behavior of the product G(rj) < p{rj) > under a variation of the observational 
scale is related to a competition between anisotropy and density contrast, and it may be of 
significance in the correct interpretation of recent cosmological data. In particular, a sufficiently 
large spatial anisotropy (even maintaining local homogeneity) may increase G{rf) < p(rj) > in 
a significant way as r) is increased. 

6.3 The choice of the initial data set 

The general picture arising from the above analysis is that we pick up an appropriate initial 
data set which, when evolved, gives rise to the real space-time. The description of this data set 
and of the resulting space-time is too detailed for being of relevance to cosmology. Intuitively, 
one would like to eliminate somehow all the unwanted (coupled) fluctuations of matter and 
space-time geometry on small scales, and thus extract the effective dynamics capturing the 
global dynamics of the original space-time. The possibility of actually implementing such an 
approach is strongly limited by the fact that we do not know a priori the structure of the 
space-time we are dealing with. But we may alternatively decide to handle the unwanted 
fluctuations at the level of data sets, since the time evolution of the initial data set for the 
Einstein equations is actually determined by the very constraints which that data set has to 
satisfy. As we have seen above, this can be done quite effectively in the framework of the 
Renormalization Group which is naturally well suited to this purpose. However, and here we 
come to the point we wish to make clear, different initial data sets giving rise to the same 
inhomogeneous and anisotropic space-time, may yield smoothed data set giving rise to different 
FLRW space-times. In other words, the renormalization procedure and the dynamics 
do not commute. The dynamics gives rise to the crossover between different FLRW space- 
times, or more generally, between different renormalized models of the same original irregular 
space-time. 

This situation is in fact not so paradoxical as it may seem. From a thermodynamical point 
of view, we have seen that one of the members of the initial data set, namely the three-metric 
g, plays the role of a temperature. Thus, by varying the metric, one can move through the 
possible pure phases of the thermodynamical system considered. In this sense, a real, locally 
inhomogeneous universe, is to be considered as akin to a generic complex thermodynamical 
system. The possible locally homogeneous cosmological models, arising from it by suitable 
choices of initial data set to smooth-out, correspond to its distinct pure phases. The resulting 
dynamics yields, in an analogy with common statistical system, a dynamical crossover between 
different pure phases. This makes accessible in cosmology too, the whole subject of critical 
phenomena with a plethora of interesting consequences. Critical phenomena are always man- 
ifested macroscopically, as phase transitions are collective phenomena in their nature. This 
aspect may turn out to be of importance for the study of structure formation and clustering in 
the universe. 
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The issue that now needs consideration concerns the existence of natural initial data sets, 
in the real lumpy universe, with respect to which the above Renormalization Group smoothing 
can be implemented. As remarked previously, the formalism is particularly well suited to a 
data set supported on slices of constant extrinsic time, (viz., spatially constant rate of volume 
expansion, k). This is so simply because the Ricci-Hamilton flow, characterizing the scale- 
dependence in the fluctuations of geometry, is most conveniently normalized to preserve volume. 
However, depending on the particular geometry we wish to smooth, different normalizations 
can be envisaged too [25], and the relevance of constant extrinsic time initial data set can be 



simply traced back to their importance in the standard analysis of the initial value problem in 
relativistic cosmology. 

A suitable slice of a frame comoving with matter is another natural choice that almost 
immediately comes about. However, in general, we cannot choose surfaces orthogonal to the 
cosmic fluid flow lines, since if rotation is present, they do not exist. A more proper choice 
would be to select the surface of constant matter density as the surfaces of constant time, 
(see [15 [ and references quoted therein). It has been argued recently that this choice is rather 



optimal. In such a slicing jLjj], and in the observable domain, the spacetime metric takes the 
form 

g ( 4 ) = -A 2 (x a )dt®dt + S 2 {t,x h )f ij {t,x h )dx i ®dxi, (158) 

where a = 1,2,3,4, i,j = 1,2,3. From an observative point of view, the function A 2 (x a ) 
is nearly constant since the acceleration of the corresponding timelike congruence is small; 
S 2 (t, x h ) is nearly independent of time since the expansion of the congruence is almost spatially 
constant, and finally also fij(t,x h ) is nearly independent of time, because the shear is nearly 
zero. These remarks suggest that the slices of constant time characterized in this way are, at 
least in our observative domain, nearly spaces of constant curvature, so that they are the most 
suitable for implementing the above smoothing procedure. 

In this connection notice that, for an initial data set associated with such a surface of 
constant matter density, some of the formulae of the previous sections are simplified. We refer 
in particular to the expression (|158| ) providing the scale dependence of G{rj) < p(rj) > which 
now reduces to 



G(v) < P(V) >= G( V o) < p(vo) > +4^ < >M &n + 0(S V 2 ). (159) 

Thus, for an initial data set associated with a surface of constant matter density, the 
product G(rj)p(yj) increases with the averaging scale as we average-out the local anisotropics 
in the geometry. And, as remarked in the previous section, Gp — Gp can be quite large, of the 
same order of magnitude of Gp, if f(x h ,t) is sufficiently spatially anisotropic. 

Finally, one may wish to choose initial data set corresponding to a frame minimizing the 
anisotropies in the CMB radiation. This is characterized by that unique four-velocity which 
eliminates the CMB dipole. This corresponds to a family of observers moving with the back- 
ground radiation. This choice has an advantage that it can be accurately determined by local 
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observations, but it is manifestly not very suitable to use in the Ricci-Hamilton formalism, as 
developed here. 



7 Concluding remarks 

According to the contents of this paper the key idea on which our whole analysis rests is that 
of RG, namely, the involved physics is that of the running (scale dependence, be it energy or 
momentum scales) of the couplings and the relevant quantities, accordingly. This philosophy 
is recognized and established in particle physics, e.g. it is well known that the fine structure 
constant a measured at low energies is different from the one measured at the LEP energy scale. 
In each case it is the presence of "fluctuations" (of any kind) that requires a scale dependent 
redefinition ("dressing") of the physical parameters which can, in turn, modify them, as well as 
the very structure of the theory, in a non-trivial way. Applications of the RG to particle physics 
have usually been in the ultra-violet limit (e.g. in QED, QCD, GUT) whereas in condensed 
matter physics they have been in the infra-red limit, in the study of critical phenomena and 
phase transitions. 

We have taken this infra-red direction in cosmology. The application of the concept of 
running of the physical quantities, motivated by RG, appears to be a new important feature 
in a cosmological setting, providing (at least) partial explanation of some controversies of 
standard cosmology, which we are going to discuss below. Let us first point out that generally 
running, also, of cosmological quantities is as such motivated by the asymptotically free higher 
derivative quantum gravity, according to which the gravitational constant is asymptotically 
free @. Taking into account this fact, of running G, one can explore its consequences in the 
standard FLRW cosmology (cf. jl9|, |5|). 

A word of caution is in order - since we do not have an ultimate theory of quantum gravity 
(QG), such approaches to cosmology are not on a rigorous basis from a theoretical point of view 
and should rather be taken as phenomenological. In principle, once we have a valid QG theory 
will one be able to directly derive a RG equations for various cosmological quantities. Although 
we have taken here a more standard view in which a split exists between the background 
(associated with infra-red effects) and renormalization of fluctuations, there may well be scales 
where such a split is not sensible at all. However, a RG capable of interpolating between 
the qualitatively different degrees of freedom in a parameter space of gravity, remains to be 
developed which may, after all, as well be possible after QG theory is within our reach. 



One of the major issues in modern cosmology is concerned with the value of the Hubble 
constant and the apparent conflict between the observed age of the universe and the predicted 
one, in the standard FLRW model, based on the recent measurements of the Hubble constant. 
Namely, the recent measurements of the Hubble constant using the Virgo cluster (distance ~ 15 
Mpc) strongly support that the Hubble constant H a has somewhat larger value h = 0.87± 0.070 
pp| . At the same time, other distance indicators yield a systematically smaller value of H , 
e.g. around 0.55 ± .08 using NGC 5253 at a distance of 4 Mpc p5|], w hile an analysis of the 
gravitational lensing of QSO 0957+561 indicates h = 0.50 ± 0.17 [41]. When one calculates 
now the age of the universe, using the larger value of H Q , in the FLRW model one runs into a 



3 /i is H a measured in units of 100 km/sec/Mpc 



49 



serious problem, as the predicted age turns out to be too small to accomodate the measured 



ages (~14-18 Gyr) of the globular clusters in our galaxy [j44|, |32f| . 

Moreover, typical inferred values of the density parameter Q Q = p /Po nt (Po is the present 
value of the total energy density of the universe and p c J lt its present critical energy density, 
defined as p^" lt = 3H% / \8irG) where G is Newton's gravitational constant) increase correspond- 
ingly with the increasing size of various structures (e.g. |29|]). These measurements can at most 
account for a fraction of 0, o which, according to the inflationary paradigm, should be equal to 
1. This in turn is one of the reasons for postulating the existence of non-baryonic Dark Matter 
(DM) which is also required to explain the structure formation. 

Various people have since then looked at possible theoretical alternatives to these DM 
scenarios, such as e.g. introducing a cosmological constant in the Einstein equations or ad hoc 
modifications of the usual theory of gravity. The important point in this respect may as well be 
the one addressed in this paper. It is usually taken for granted that, on large scales, the universe 
is described by the FLRW solution. There is no alternative really since we do not know any 
solution of Einstein's equations capable of describing a clumpy universe. Nevertheless, even in 
the absence of explicit fine-grained models, we would like to know how in principle, and when, 
one could extract a background model from an inhomogeneous one, such that: (i) they both 
obey, "approximately", the Einstein equations despite the averaging or smoothing involved, 
and (ii) observational determinations of cosmological parameters (H , o ,...) correspond in 
a sensible way to that mathematical averaging procedure. Thus an issue of importance for 
cosmology [37], is the question on what scale is the FLRW model supposed to describe the 



universe? Likewise, what averaging scale are we referring to when we give the value of Q Q , 
whose definition necessarily refers to an idealised, i.e. smoothed background model. 

There has been recently an increased effort in this direction with some interesting results, as 
e.g. that the coarse-graining effects could be non-negligible in the context of affecting the age 
of the universe. For example, considered a model, of locally open (underdense) universe 
embedded in the spatially flat universe, in which the expansion rate in our local universe is 
larger than the global average. Similar model was considered in |38|] where a local void in the 
global FLRW model was studied and the inhomogeneity described by the Lemaitre-Tolman- 
Bondi solution. The results indicate that if we happened to live in such a void, but insisted on 
interpreting the observations by the FLRW model, the Hubble constant measurements could 
give results depending on the separation of the source and the observer, providing a possible 
explanation for the wide range of their reported values and capable of resolving the age-of- 
universe problem. 



On the other hand [19] studied the QG effects at cosmological scales (the phenomenon in 
question is exactly that of quantum coherence, known from a laboratory to happen on macro- 
scopic scales of the order of cm), assuming asymptotic freedom of the gravitational constant 
and incorporating running G, according to the appropriate RG equations, into FLRW model. 
Such G takes the value of Newton's constant Gn at short distances but then slowly rises as 



distance increases. However, as mentioned in [19], the RG equations used there might not be 



applicable in the infra-red regime studied, but these concerns were put aside having in mind 
the absence of any other available beta functions for QG in the infra-red regime. 

One can also approach the averaging problem, modifying the FLRW metric (or equivalently 
Cosmological Principle) and the Einstein equations, by an introduction of a generalized scale 
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factor which depends both on t and the scale r [28]. This introduces the scale dependence of 
(Gp) and running of other cosmological quantities such as H Q , fl Q , and the age of the universe 
t Q , as functions of distance scales. 



Let us notice that this picture of running of cosmological quantities comes about natu- 
rally in our approach which is physically motivated by RG, namely, due to increasing of the 
gravitational constant with scale, (and possibly increasing amount of DM), as discussed in the 
previous section, £1 has effectively increased too. Moreover, since the scale factor is governed 
by the scale dependent (Gp), it now seems to depend on the scale r as well, i.e. it increases 
at any fixed time with the increasing distance. Consequently, the value of H Q is not the same 
everywhere in the observable universe and depends on a scale where it is measured. Moreover, 
since H Q ~ l/t Q the universe becomes older when its age is estimated on a smaller scale. This 
is not to be taken as implying that the age depends on where one calculates (every observer 
using the same scale r at some time t\ will obtain the same age). The key point to emphasize is 
that, having in mind the RG arguments and interpretation, a direct comparison of cosmological 
quantities makes sense only when they are measured (or calculated) with respect to the same 
scale, since the same quantity can take different values at different scales. 

Notice that independently, also Quantum Cosmology advocates, though in a different con- 
text of bubble universes, a possibility that we may live in a universe in which the value of 
Hubble constant and the measured density are different in different places and in our local 
neighborhood Q Q may well be less than 1 [j34f . 

This makes our proposal even more promising. 
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Figure Captions 



Fig. 1 A portion of a minimal geodesic balls covering. The dotted disks are the |-balls 
which pack a given region. The larger, undotted disks, represent the e-balls which provide the 
covering. 

Fig. 2 We can approximate the average over E by averaging over euclidean balls whose 
Lebesgue measure has been locally weighted through Puiseux' formula. In this drawing, the 
^-geodesic balls are represented by curved disks on E, while the corresponding euclidean balls 
are correctly depicted as three-dimensional. 

Fig. 3 The local average of a function / over a geodesic ball B(x, e G ) feels the underlying 
curvature of the manifold through Puiseux' formula. In particular, in passing from B(x,e ) 
to the larger ball B(x,e + rj) we get correction terms which depend on the fluctuations in 
curvature. 



Fig. 4 The intricate but symmetrical intersection pattern of an array of |-balls in the 
2D-plane. Here, the |-balls are dotted while the overlapping e-balls, providing the covering, 
are the dashed circles. The solid circles and the arcs of circles describe the intersection pattern 
of a |e-ball. On a curved manifold of bounded geometry, the pattern is more complicated and 
not symmetrical at all. Nonetheless, one may easily obtain a recursive definition of blocking 
by exploiting a partition of unity argument, explained in the text. 

Fig. 5 The intersection pattern of a 2e-ball B(xj,2e) with five e-balls B(xh,e). The 2e- 
ball is represented by a dotted disk, while the e-balls are represented by solid circles. The 
white disks represent the packing of B{xj, 2e) with |-balls. The standard representation of the 
integral over B(xj,2e), through the partition of unity subordinated to the B(xh, e)-covering, 
yields the recursive relation ip m+ i(j) = J^h^mih)- 

Fig. 6 Even if the function f is constant over E, its average values over the geodesic 
balls, B(xi, e) and B(x2, e), are generally uncorrelated owing to the random fluctuations in the 
geometry of the balls. 

Fig. 7 The matter distribution p, J, the lapse function a, and the shift vector field a 
associated with the istantaneous observers on E. Here, n denotes the four-dimensional unit, 
future-pointing, normal to E. The vector field an + a connects the corresponding points in the 
two infinitesimally close slices E and E^. 
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Fig. 8 Given a probability law, according to which matter is distributed at a given length 
scale (say, the planetary scale), we can get the corresponding probability distribution obtained 
by averaging the matter variables over regions of ever increasing scales. The resulting effective 
Hamiltonians are defined up to additive constants which affect the renormalized mass den- 
sity and the renormalized momentum density. Such indeterminacy can be naturally removed 
by enforcing the Hamiltonian and divergence constraints at each step of the renormalization 
procedure. 



Fig. 9 Upon enlarging the ball, from e to e 7 e, we alter the average value of the matter 
variables p and J, (even if p and J do not vary appreciably), since there can be strong curvature 
fluctuations in the annulus B(x,e 1 e)\B{x,e). 



Fig. 10 If we enlarge the ball from e to e v e, while deforming the metric according 
to the Ricci-Hamilton flow (52), (for a parameter time rj = 7), then the average value re- 
mains as constant as it can be since (|52]) smoothes the curvature fluctuations in the annulus 

B(x,e 1 e)\B{x,e). 



Fig. 11 If we blow up the metric of two three-spheres and S%\, (here represented as 
two-spheres), around the points y\ and 1/2 and then join the resulting manifolds, we obtain the 
connected sum S^AS^y According to this procedure, the neck S 2 x IR 1 inherits a cylindrical 
metric up to exponentially small correction terms. 



Fig. 12 The Ricci-Hamilton flow of S^AS^ can be- expressed in terms of the ^-evolution 
of the inclusion maps X(»)j ct = 1,2, 3, and of the known evolution of the locally homogeneous 
metrics on S9 i AB(y i ) and S 2 x IR 1 . 



Fig. 13 Different deformations of the initial three-manifold £ may have quite different 
fates. If the waist of £ is increased enough by the deformation, (i.e., if it is rounded), the 
Ricci-Hamilton flow will evolve £ toward the round three-sphere. If the waist is shrinked 
enough by the initial deformation, the Ricci-Hamilton flow will evolve S toward S%\ U S%\ . 
The three-manifolds attracted towards Sf^ ]J Sf 2 \ define a critical surface. 



Fig. 14 Two neighbooring three-manifolds, (£, g) and (£, g + SK), with the same volume, 
evolving under the Ricci-Hamilton flow toward isometric round three-spheres (S 3 ,g can ) and 
(S 3 ,g* an ). In general, g can and g* can differ by an infinitesimal diffeomorphism </> generated by a 
vector field v, i.e. g* an = g can + 5 [Lv9 can — ^ffcan^jf*^ . By exploiting the freedom in choosing 
the shift vector field a, we can get rid of (p. 
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Fig. 15 The time evolution of an inhomogeneous and anisotropic initial data set and 
its Ricci-Hamilton renormalized counterpart. Notice, that (M 4 ,g( 4 )) depends in a sensible 
way on the particular initial data set (£, g, K, p, J) which is renormalized according to the 
Ricci-Hamilton flow. Different initial data sets, corresponding to the same inhomogeneous 
and anisotropic space-time (M 4 ,^ 4 )), may generate distinct renormalized model space-times 
(M 4 ,g( 4 )). This seemingly paradoxical situation has a natural explanation in terms of the RG 
approach. 
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